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ABSTRACT 


A package of general and efficient computer code 
for solving optimization problems is presented in this 
thesis. Powerful unconstrained optimization methods 
based on the ideas of quadratic convergence and conjugate 
directions are used. Constrained optimization problems 
are solved, using the penalty functions and feasible 
directions method. 

These programs are tested for convergence and 
efficiency by using typical test functions previously 
suggested in the open literature. Results of a compara- 
tive study of constrained optimization methods indicate 
that sequential penalty function methods, with a gradient- 
based search procedure, using finite difference gradient 
computations, perform well on most problems. 

The capability of the package is well demons- 
trated by the design of a direct transfer air-to-water 
compact heat exchanger employing flat-finned tube surface, 
the objective being minimum pumping power and the cons- 
traints imposed on volume, dimensions of the exchanger 
and heat transfer. The £ -NTU method is used for evaluat- 
ing the heat exchanger performance. Polynomial expressions 



are found to approximate the experimental data for the 
fluid transport properties and the friction and colburn 
factors for the heat transfer surface. This optimization 
study of the compact heat exchanger provides a design 
which has got much superior performance. 
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CHAPTER 1 


INTRODUCTION 


In modern technology, it is becoming increasingly 
important for engineering designs of all types to offer 
the best possible performance- To achieve this task, the 
analyst must identify the variables which influence the 
behaviour of the system, judiciously select quantitative 
measures that appropriately describe the system perfor- 
mance and then make a detailed parametric study. Fre- 
quently, it is found that some of the desired objectives 
have conflicting trends. Also, in many applications it 
is not possible to keep the design variables constant as 
they may have a tendency to drift with time or fluctuate 
between some limits. For a reliable engineering design 
it is not adequate to just determine the best operating 
conditions. It is equally important to know how sensitive 
the optimum performance is/to changes in the system para- 
meters. Another important consideration that must be 
kept in mind while making a design study is that the 
design requirements may change in future, and it would 
be very useful if a suitable alternate design can be 
predicted from the results of the current optimum. This 
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would eliminate the need for making another extensive 
parametric study. This necessitates the introduction of 
a design tool that will permit the engineer to experiment 
with design variables and different performance criteria, 
to search for the best operating conditions and to identi- 
fy the influence of the system parameters as the selected 
performance indices. It is also desirable to include 
some aids to show how the design specifications constrain 
the system performance. 

Computer-aided design procedures are playing an 
important role in the design practice, making it possible 
to analyze a wide variety of complex problems. Recent 
advances in computer-graphics, real time manipulation 
capabilities along with the developments in optimization 
theory and interactive computation have made it possible 
to tap the design tools'potential . The present work is 
concerned with the development of a single package which 
includes several optimization methods. The package has 
been used for optimum design of a compact heat exchanger. 

1.1 HEAT EXCHANGER OPTIMIZATION 

Heat exchangers are designed for many different 
applications and hence may involve many different perfor- 
mance criteria. These criteria may be: 

minimum initial cost, minimum initial and 
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operating costs, minimum weight or material, 
minimum volums or heat transfer surface area, 
and so on. 

When a single performance measure has been defined quanti- 
tatively and is to be minimized or maximized, it is called 
an "objective function" in a design optimization. A parti- 
cular design may be subjected to certain " CONSTRAINTS" 
such as desired heat transfer rate, allowable pressure 
drop, limitations on height, width and/or length of the 
exchanger. A number of different surfaces could be 
incorporated in a soecific design problem and there are 
many geometrical parameters that could be varied for each 
surface geometry. For an extended surface exchanger, the 
geometrical variables associated with a fin are the fin 
pitch, fin height, fin thickness, type of fin, and other 
variables associated with each fin-type. In addition, 
operating flows and temperatures could also be changed. 

Thus a large number of " design variables" are associated 
with a heat exchanger design. The question arises as to 
how one can effectively adjust these design Vc a-iauies 
within imposed constraints and come with a design hav- 
ing optimum objective function. This is what we mean by 
the most " efficient" design. 

A complete mathematical optimization of heat 
exchanger design is neither practical nor possible. 
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Many engineering judgements based on experience are in- 
volved in different stages of the design. However, once 
the general configuration and surfaces are selected for 
a specific application, an optimized heat exchanger design 
may be arrived at if the objective function and constraints 
can be expressed mathematically, and if all of the varia- 
bles are automatically and systematically changed on some 
statistical or mathematical basis. 

A large number of optimization (search) techniques 
are available in literature. The difficulty with the pre- 
sent state of optimization art is that while any single 
search technique may work well on many problems, no single 
technique is able to solve every problem. A package of 
optimization programs is a convenient approach to over- 
coming this difficulty. Well developed, the programs can 
be used interchangeably and a new search can be called 
in, where a previous one has quit thereby continuing or 
assuring progress toward optimum. In addition, when a 
package is used, various parts of the search procedure 
common to several optimization routines need be programmed 
only a single time. 

1.2 OPTIMIZATION 

Optimization techniques offer systematic approach 
for seeking the best operating conditions. A description 



of many of the more popular optimization algorithms is 
presented in [l, 2, 3 ]. Almost all design problems require 
either the maximization or the minimization of some para- 
meter. This parameter is usually called the design objec- 
tive function. For example the problem may call for a 
heat exchanger with a minimum volume. The expression for 
volume would be the design objective function. This is 
done by selecting values for certain variables. These 
variables are known as design or decision variables. For 
the design to be acceptable/ it must satisfy certain cons- 
traints. For example, an air heater must be designed so 
that it will fit into a given space. 

The general nonlinear constrained optimization 
problem can be written mathematically as: 


Minimize 

F ( 

x ) 



(1) 

Subject to: 

g i ( 

X ) < 

0 

i = 1 , 2 , . . ,m 

(2) 


V 

X ) = 

0 

k = 1 , 2 , . . ,p 

(3) 


X* < 

q - 

- *q - 

x u 

q 

q = 1 , 2 , . . . , n 

(4) 


Where the vector, X , is the vector of " n " design vari- 
ables, F is the objective function, g^ is a set of " m " 
in~equality constraints, 3^ is a set of " p" equality 
constraints and the lower and upper bounds , X^ are 
the side constraints imposed on the design variables. 

Side constraints could also be included in (2), but for 
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better efficiency they are treated separately by some 
software experts. Equations (2) and (3) may be linear 
or nonlinear functions of the design variables. They may 
be explicit or implicit functions of X , but must have 
continuous first derivatives. The basic solution proce- 
dure is to, systematically, search the design space by 
sequentially varying the design parameters along search 
directions that reduce the objective function. For 
example, starting at a point defined by the existing 
values for the design variables, a new point X is 
determined as: 



X + a s 

q q 


F< V > 4 F( V 

Where S is the search direction, and or. is the step size 
along that direction. The essence of nonlinear programm- 
ing methods lies in determining the oc which minimizes 
F ( X by using a one-dimensional search technique. A 

new search direction is computed at X^ + ^. A minimum is 
obtained if it satisfies the accuracy criterion. There 
is no evidence that the resulting minimum is the global 
optimum, since relative minima may be present. The search 
procedure must be repeated using different starting points 
if this is a possibility. For problems based on real 
processes, luckily, the objective function is usually a 
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well behaved function with a single extremum. Therefore/ 
for most practical purposes / the use of standard numerical 
procedures that provide a local solution to the optimiza- 
tion problem is not a great disadvantage. 

1.3 PREVIOUS WORK 

Compact heat exchanger design is a complex task 
requiring the examination and optimization of a wide 
variety of heat transfer surfaces. Smith [ 4 ] has listed 
some typical advantages of direct cooling with air as 
compared to cooling with water in a shell-and-tube 
exchanger. Studies have shown that a poor choice of 
either the heat transfer surfaces or design parameters 
can more than double the costs chargeable to a heat 
exchanger. 

The method of Bergles et al. [ 5 1 is primarily 
meant for comparing the performance of heat exchanger 
surfaces for some specified criteria by adjusting only 
two exchanger variables. Their method does not include 
any minimization technique but results show that a great 
improvement in heat exchanger performance can be made by 
proper selection of design parameters. 

Pax and Mills [e] developed a method, using 
Lagrange multipliers to optimize heat exchanger design 
under specified constraints. This technicrue requires 
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the objective function and the constraints expressed in 
explicit equation form and differentiable throughout the 
range of interest/ the total number of constraints be 
less than the total number of variables, and all cons- 
traints be equality constraints. Obviously, the use of 
this method is restricted to a limited number of problems. 
Briggs and Evans [?] recommended a " logic search method" 
in which an experienced engineer makes use of selecting 
design variables based on his apriori knowledge to obtain 
an optimum design. This method is less scientific and 
can be tried only if the computer time and storage capa- 
city are freely ensured. Wilde [3] has surveyed various 
optimization techniques applicable to heat exchanger 
design. He focussed on the geometric programming tech- 
nicrue and other search methods which require explicit 
expressions for the constraints and the objective func- 
tion. The geometric programming technique also requires 
that all of the constraints and the objective function 
must be expressed in a particular power law equation form. 
These methods have limited applications for a general heat 
exchanger optimization problem. 

Palen et al. [9] proposed the Complex Method 
[l], for the heat exchanger optimization problem. They 
found a minimum shell -and- tube exchanger by varying six 
geometrical parameters. The Complex Method requires 
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several feasible starting designs before optimization can 
be performed. Johnson et al. [lo] coupled an existing 
shell-and-tube condensor design code with a constrained 
function minimization code to produce an automated marine 
condensor design program of vastly different complexity. 

Pontein and Wassink [ll] utilized the Simplex 
Method [l] and steepest descent method for optimizing a 
shell-and-tube exchanger. Shah et al . [l2, 1?] have 
initiated work on package approach lines to the optimiza- 
tion problems in compact heat exchangers. The package 
includes the numerical nonlinear programming techniques. 
Hedderich et al . [ 1 4 3 have also developed a massive code 
using the method of feasible directions and the augmented 
Lagrange multiplier method. It has been used for the 
design of an air-to-water finned- tube heat exchanger. 

It is experienced that/ although, there are many 
methods that have been presented for heat-exchanger opti- 
mization, each of the methods has its ox^m limitations; 
none is completely general. Among the design procedures 

cited above. those which are applicable . to cross flow, 

) 7 

air cooled compact heat exchangers are limited to the 
120 individual surfaces found in the open literature C 13 3 
for the calculation of the air side heat- transfer co- 
efficient and friction factor. Therefore, the designer 
is faced with choosing an optimum surface from a number 
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of individual optimal designs calculated from one of the 
above methods. 

1.4 PRESENT WORK 

The work embodied in this thesis is an attempt 
in two directions: 

(i) to develop a general purpose computer aided 

design procedure that assists the engineer in 
establishing a best design. This makes avail- 
able to the user, numerous nonlinear programming 
technigues. The search parameters can be changed 
and the user can switch from one search technique 
to another. This flexibility, provided to the 
user, gives a feel for the various technique's 
applicability and to select the most suitable 
one to the problem in hand. The development of 
this computer package is described in the next 
chapter, 

(ii) to use the above package for optimizing the 

design of a compact heat exchanger for some 

specific objective and constraints. A direct 

transfer, air-water finned- tube heat exchanger 
1 *) 

for a gas turbine plant has been optimized. The 
objective is to minimize pumping power require- 
ments of the heat exchanger which has a volume 
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less than a specifier) value and for which the 
heat transfer rate is prescribed. The flow 
rates of the fluids and the linear dimensions 
of the exchanger comprise the design variables 
of the problem. Fixed bounds on these variables 
are also imposed as side constraints. Cons- 
traints take into a account the piemsthetic 
value to get a realistic configuration for the 
exchanger. 

The package implements the following multi- 
dimensional minimization schemes: univariate, steepest 
descent, conjugate directions, conjugate gradient, 
variable metric method [1,2,3]. These methods use either 
a quadratic interpolation procedure or golden section 
technique [l, 2 ,?] for 1 -D minimization. Constrained 
problems are solved using either penalty function or 
feasible directions [l,?] method. To check the relia- 
bility and flexibility of the proposed package, a number 
of problems available in literature have been solved and 
results compared and listed in Appendix-A. These problems 
differ fundamentally in- tbe manner in which the objective 
functions and th« constraints- ax-e compute*! in tire- func- 
tion subroutine. These quantities are usually programmed 
±p. one user written function subroutine as follows: 
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1. Closed form algebraic expressions 

2. Mostly algebra but requiring some iterations 

The programs are simple to use, with the user having to 
provide a function subroutine and to specify a minimum 
number of variables and flags. The ensuing chapter 
deals with the package structure that has been developed. 



CHAPTER 2 


PACKAGE STRUCTURE 

2 . 1 INTRODUCTION 

A set of general programs have been developed 
which make available to the user, reliable mathematical 
programming routines suitable for optimizing a large 
class of problems. Various numerical techniques which 
have been utilized to structure this optimization package 
are briefly discussed in this chapter and summarized in 
Pig. 1. Flow diagrams for solving unconstrained and 
constrained optimization problems are given in Fig. 2. 
Complete details regarding the algorithms used, may be 
obtained from [ 1 , 2 , 3 ]. Appendix-A gives the details for 
using the programs with examples to illustrate the proce- 
dure and program-listings have b<=en given in Appendix-B. 

An interesting feature of the package is it's 
flexibility. New search methods can be easily introduced 
by making minor changes. This is made possible by having • 
maximum interaction between the subroutines and dividing 
the procedure into the following four basic operations; 

1. Selection of a search direction 

2 . Estimation of the range for minimum along 
the search direction 
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Fig 2(b). Flow chart for Constrained Optimization 
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?. One-dimensional minimization to locate the 
minimum along the search direction 

4. Checking for convergence 

The multidimensional methods differ primarily in the 
manner by which the search direction is determined. The 
effectiveness of the methods for computing a given search 
direction depends on the accuracy with which the minimum 
along previous search directions has been obtained. Thus, 
the one-dimensional minimization influences the stability 
and convergence of the multi-dimensional methods. 

2.2 RANGE-FINDER: SUBROUTINE RANGE 

Nonlinear programming methods essentially con- 
sist of a series of unidimensional minimizations. Numeri- 
cal solution procedures for locating a minimum in a search 
direction require a prior knowledge of the limits between 
which the minimum exists. These limit points are deter- 
mined by the range-finder by sequentially increasing the 
step-size along the search direction until three consecu- 
tive points are obtained, such that 

F( *q + “l V > F( *q + a 2 V 

(5) 

F< X q + a 3 y > F( x q + a 2 y 

where a. , or , a_ define three consecutive points 
±23 

along the search direction and X g is the design vector 
the S is search direction. The Equation (5) guarantees 

q 
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that a minimum exists between a and a^. 

The inputs to RANG E are X g , S g , initial step- 
size cl, maximum increment in stepsize and the number of 
iterations, permitted to determine a , a 2 , and a^. It 
is desirable to obtain a small value for (a^ - ) in a 

few iterations. If a step along S reduces F, a is 
doubled for the next iteration, and the values for cc^ and 
(*2 updated. The increment in a is not allowed to exceed 
the maximum specified value. When more iterations are 
required to estimate (oc^ - ct^) than that specified, RANGE 
assumes the last value of a to be the location of the 
minimum along S . While seeking a 9 , RANGE has the pro- 
vision to reverse the search direction and reduce ct, if 
an initial step along S increases F. 

Since RANGE is used also by constrained optimi- 
zation methods, some additional features have to be intro- 
duced to take care of constraints. When using the method 
of feasible directions, it is necessary to locate cons- 
traints that are violated. If the point violates 

constraint "g" , RANGE repetitively uses a linear inter- 
polation scheme to find cl' , where 


_ % g ( W - Vn gt V 

g< x p+1 ) -g ( x p ) 


P - X q + a p S q 
Vl - + °P+1 \ 


( 6 ) 
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where g ( X p ) , g ( X p+ ^ -* are constraint function 
such that g ( X + or.' S ) <. e ~ the accuracy needed. 

ti M. 

In using the interior penalty function method, 
the search is restricted to the feasible region. When 
^ p+1 violates " g ", RANGE progressively reduces the step 
size to find <x' , which is now expressed as: 


a ' = a + 

P 


a „ - a 

P+1 p 


(7) 


such that g ( X + or/ S ) < 0 

q q — 


2 . 3 ONE-DIMENSIONAL SEARCH METHODS 


The range-finder determines the interval which 
contains the minimum along a search direction. An effi- 
cient unidimensional search is next, needed to accurately 
locate the minimum in a few iterations. The 1— D search 
problems are of the form: 

minimize 0(x) 

where x is a single variable whose optimum value is to 
be determined, so that the function 0 is at its minimum. 

The program package gives the user, the option 
of using either the quadratic interpolation scheme (SUB- 
ROUTINE QUAD) or the golden section method (SUBROUTINE 
GOLD) [l3]. The former performs a series of iterations 
approximating 0(x) as a quadratic function; the lsrtter 
isolates the minimum in regions of successively decreasing 
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size. The number of iterations needed for convergence 
to the minimum, will depend on the accuracy needed. Com- 
pared to quadratic interpolation scheme, the golden sec- 
tion method converges slowly to the minimum. 

2.4 MULTIDIMENSIONAL UNCONSTRAINED SEARCH METHODS 

Unconstrained searches, i.e., to just minimize 
F(X) without any constraints, can be performed by a sequence 
of one-dimensional minimizations in appropriate directions. 
The most efficient search procedures generate directions 
that are conjugate. 

In direct search strategies, the search directions 
are established apriori but may be modified during the 
course of a search. In the univariate method (SUBROUTINE 
UNIV ) , the design space is searched along the coordinate 
directions, every time adjusting one variable and keeping 
others constant. A cycle is completed when all the coor- 
dinate directions are searched and repeated, if needed, 
to satisfy the accuracy criterion. Powell's method 
(SUBROUTINE CONDIR) begins with the univariate approach 
but establishes " conjugate" directions by moving along 
vectors connecting points in the design space, resulting 
from previous 1-D minimizations. Conjugate directions 
[l,2,3] lead to speeded convergence of an optimization 
on most functions and, inf act, guarantee convergence on 
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an n-dimensional quadratic function after n, 1-D minimi- 
zations in mutually conjugate directions. 

Gradient based search techniques utilize gradient 
information in establishing directions for 1-D minimiza- 
tions. The objective function has a maximum rate of 
change along this direction. Thus, moving in the negative 
gradient direction would seem to reduce the performance 
index most (SUBROUTINE STEEP). The conjugate gradient 
(SUBROUTINE CONGRA) and the Davidon-Fletcher -Powell 
(SUBROUTINE DFPM ) methods [l,2,3] utilize the past and 
current gradient values to generate conjugate directions 
and speed the convergence. The Davidon-Fletcher -Powell 
approach is generally considered to be one of the most 
effective methods for unconstrained optimization. Diffi- 
culties seldom arise except on very badly distorted or 
ecentric functions. When the search direction fails to 
reduce the function, the search procedure is restarted 
after adjusting certain parameters accordingly. 

2.5 MULTI -DIMENSIONAL CONSTRAINED SEARCH METHODS 

The constraints in an optimization problem can 
be enforced either by using substitutions [ 1 , 3 ] to eli- 
minate them or by modifying the objective function [ 1 , 2 , 3 ] 
with penalty terms. The penalty terms increase the 
objective function when a constraint is violated. 
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Substitutions and penalty terms transform the constrained 
problem into a sequence of unconstrained problems. Alter- 
nate procedures for obtaining the constrained minimum are 
based on finding search directions [ 1 , 3 ] that are directed 
towards the feasible region and away from the "constraints 
just violated. The search for the optimum is confined to 
the feasible region. 

2.5.1 Interior Penalty Function Method 
SUBROUTINE IPENAL 

The basic idea underlying this method is to solve 
repetitively a sequence of unconstrained problems whose 
solutions, in the limit, approach the minimum of the 
constrained optimization problem. The objective function 
is modified as 

m 1 '1 9 

0(X,R) = F(X) -RE — — - + — -E [l . (X)J 

j=l g j(X) v/R k=l 

..... (S) 

where 0 is the modified function called PENALTY, R is 
penalty factor, " m" is the number of inequality cons- 
traints, and " 1" is the number of equality constraints. 
The penalty term with inequality constraints is small 
for points away from the constraints but blows up as the 
constraints are approached. As R is decreased, the 
penalty term, with equality constraints 1^ , blows up 
and it forces L. to be satisfied to some degree. When 

K 
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g . and 1^. are satisfied s 

0(X*,R) ->* F(5P) as R O , 

X is the optimum solution. 

An extrapolation scheme suggested by Fiacco and 
McCormick [l,2,3] is used to predict a new starting point 
and the value of the constrained optimum. The whole 
search should be strictly within the feasible region, 
because a mirror image exists in the infeasible region. 
Thus, precaution is taken in RANGE so as not to overstep 
into the infeasible region. An appealing feature of this 
method is that it produces an improving sequence of accep- 
table designs; hence allowing the designer to compare sub- 
optimal designs to the optimum one. 

2*5.2 Exterior Penalty Function Method 
SUBROUTINE EXPEN 

A frequently used exterior penalty function which 
accommodates inequality and equality constraints is of 
the form 

m o 1 ^ 

0(X,R) = F(X) + R S < g .(X) > + R 2 [ !,(£}] 

j=l J k=l K 

(9) 

where < gC^> = 0 for gjx) <_ 0 

= gjtx)for. gjC>)> 0 

0{X,R) is minimized for increasing values of R. As R 
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is increased, the terms involving g. and 1, become bigger. 

J k 

The X vector is forced by the penalty terms to satisfy 

the constraints to some degree. As long as g . and 1. are 

J K 

satisfied, as X — > X* , the value of the penalty becomes 
negligible: 

(2KX*, R) — ^ F (X*) 

X is the constrained optimum solution. This method 
converges to the constrained optimum from outside the 
feasible region for a sequence of unconstrained minimiza- 
tions with increasing R. 

Selection of proper initial values of R depends 
on the design problem and one's judgement. Fox [l] , 

Rao [ 3 ] discussed some rules for picking an initial value 
of R. A source of trouble in the penalty function methods 
lies in -the relative magnitudes of the constraints. If 
g^=1000 g£ , g^ changes more rapidly than g 2 ^hence over- 
powers it over the infeasible region. 0 is insensitive 
to g 2 f and the optimum may even violate prob- 

lem is overcome by scaling the constraints. 

2 . 5-3 Feasible Direction s Method 
SUBROUTINE FEAS 

This is based on Zoutendi jk's procedure and 
begins at a point within the constraints, using an un- 
constrained optimization procedure until a constraint 
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is encountered. Then, a usable feasible direction S is 
defined, which reduces F and does not violate the cons- 
traints, by solving the following linear programming 
problem: 

Max IS subject to S T V F + P <_ 0 
S , e S T V g + 8 e < 0 

l s il 1 1 

where 3 is a slack variable, T indicates transpose, J is 
the set of active constraints and the s^ are the compo- 
nents of S. The 0j are the " push-off " factors and, as 
suggested by Fox, are set to unity for nonlinear cons- 
traints and zero for linear constraints. 

2.5.4 Convergence Rule 
SUBROUTINE CONVRG 

For the multi-dimensional search techniques an 
optional convergence rule is provided. Besides satisfy- 
ing the accuracy criterion, this option perturbs the 
design variables and uses a scheme similar to Powell's 
[l]. It works as follows: 

1. Apply the search procedure till the accuracy 
is achieved. Call the resultant point A. 

2. Perturb the design variables by lj 


| do) 
j euj 

i = 1,2 , . . ,n 
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3. Apply the normal search until the accuracy is 
achieved. Call the resultant point B. 

4. Terminate if F(a) and F(B) differ by the accuracy 
desired, or else search for minimum along the 
line joining A and B and starting from a point 
having lower performance index (say A); call it 

C and neglect B. 

5. Terminate if F(C) and F(a) satisfy the accuracy 
criterion, or else reoeat step 1 using C as 
starting point. 

The above scheme is less likely to stop prematurely, but 
it may prove to be time consuming. 

Chapters 3 and 4 describe a compact heat exchanger 
and it's design methodology to help understanding how the 
optimization package, developed in this chapter, may be 
used for optimizing the design of a compact heat exchanger. 



CHAPTER 3 


COMPACT HEAT EXCHANGER 

3.1 INTRODUCTION 

A heat exchanger is a device which provides for 
transfer of kKfe'xxai thermal energy between two or more 
fluids at differing temperatures. Heat transfer between 
fluids takes place through a separating wall. Since the 
fluids are separated by a heat transfer surface, they 
do not mix. Common examples of such heat exchangers are 
the shell-tube-ex changers, automobile radiators, condensers., 
evaporators, air preheaters. A heat exchanger consists 
of the active heating elements such as a core or a matrix 
containing the heat transfer surface, and passive fluid 
distribution elements such as headers, manifolds, tanks, 
inlet and outlet nozzles or pipes, or seals. Usually, 
there are no moving parts in a heat exchanger; however, 
there are exceptions such as a rotary regenerative ex- 
changer, in which the matrix is mechanically driven to 
rotate at some design speed. 

The heat transfer surface is the surface of the 
exchanger core which is in direct contact with the fluids 
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and through which heat is transferred by conduction. That 

portion of the surface which also separates the fluids is 

referred to as " primary or direct surface". The design 

thermal effectiveness of exchangers employing such surfaces 

is usually 60# and below and the heat transfer surface area 

2 3 

density is usually less than 300 m /m . In many applica- 
tions, a much higher, upto 93 % , exchanger effectiveness 
is essential, and the box volume and mass are limited 
so that a much more compact surface is mandated. Usually, 
either a gas or a liquid having a low heat transfer co- 
efficient is the fluid on one or both sides. This results 
in a large heat transfer surface area requirements. For 
low density fluids (gases), pressure drop constraints, 
tend to require a large flow area. So a question arises, 
how can we increase both the surface area and flow area 
together in a reasonably shaped configuration. Appendages 
or fins on the primary surfaces increase the surface area 
density. Flow area is increased by the use of thin gauge 
material and sizing the core property. The heat transfer 
coefficient on the extended surfaces may be higher or 
lower than that on the unfinned surfaces. For example, 
the interrupted (strip, louver, etc.) fins provide both 
an increased area and increased heat transfer coefficient, 
while the internal fins in a tube, meant primarily for 
structural strength and flow mixing purposes, may result 
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into a slight reduction in the heat transfer coefficient 
depending on fin spacing. 


Loosely defined, a compact heat exchanger is one 

which incorporates a heat transfer surface having a high 

area density". That is, a high ratio of heat transfer 

surface to volume. Some what arbitrarily, it is specified 

that a compact surface has an area density 3 greater than 
2 3 

700 m /m . A spectrum of surface area densities of heat 
exchanger surfaces is shown in Fig. 3. On the bottom of 
the figure, two scales are shown: the hydraulic diameter 
in mm and equivalent heat transfer surface area den- 
sity 3 (m /m ). Different exchanger - surfaces are shown 
in rectangles. The short vertical sides of a rectangle, 
when projected on the 3 (or ) scale, indicate the range 
of surface area density (or hydraulic diameter) for the 
particular surface in question. Interesting details, for 
heat exchanger classification have been given [l6] with 
adequate figures. The motivation for using compact sur- 
faces is to gain specified heat exchanger performance, 
q/^t m , within acceptably low mass and box volume cons- 
traints. As 


g_ 

At 


m 


u 3 v , 


where q is the heat transfer rate and At 
A m 

temperature difference. 


the mean 
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Clearly, a high P minimizes the exchanger volume V. 
Moreover, compact surfaces generally result in a higher 
overall conductance U. This also amounts to having a 
smaller volume. As compact surfaces can achieve struc- 
tural stability and strength with a thinner section, the 
gain in lower exchanger mass is even more pronounced than 
the gain in lower volume. Various techniques employed 
to make heat transfer surfaces compact are: fins between 
plates, finned circular tubes, or densely packed conti- 
nuous or interrupted cylindrical flow passages of various 
shapes . 

A heat exchanger of any structural construction 
is considered as compact if it employs a compact surface 
on either one or more sides of a two-fluid or multi-fluid 
heat exchanger. The convective heat transfer coefficient 

for gaseous fluids is generally one or two orders of 

Sfose c V 

magnitude lower than^water, oil, and other fluids. Thus, 
to reduce the size and weight of a gas-to-liquid heat 
exchanger, the heat transfer surface on the gas side 
needs to be much more compact as compared to a simple 
heat exchanger having only circular tube. Hence, for a 
somewhat " balanced" design, a compact surface is employed 
on the gas side. Thus major applications of compact heat 
exchangers are gas-to— gas, gas— to— liquid, and gas— to- 
condensing or evaporating fluid heat exchangers. 
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3.2 CONSTRUCTION TYPES AND SURFACE GEOMETRIES 

The basic construction types employed in the 
design of a compact heat exchanger are extended surface 
heat exchangers employing fins on one or more sides, 
regenerators employing small hydraulic diameter surfaces, 
and tubular exchangers employing small diameter tubes. 

Two most common types of extended surface 
exchangers are the plate-fin and tube-fin types. In a 
plate-fin exchanger, fins or spacers are sandwiched 
between parallel plates, referred to as parting plates 
or parting sheets or formed tubes. Fins are attached 
to the plates by brazing, soldering, gluing, welding, 
mechanical fit, or by extrusion. While the plates 
separate the two fluid streams, the fins form the indi- 
vidual flow passages. Alternate fluid passages are 
connected in parallel by suitable headers to form two 
or more sides of the exchanger. The plate-fins are 
categorized as 

1. Plain (uncut surfaces) and straight fins 

2. Plain but wavy fins 

3. Interrupted fins 

Plain triangular and rectangular fins obviously belong to 
first category. Strip, louver, perforated and pin fins 
comprise the last category. The heat transfer coefficient 
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and friction factor in the developing flow region are 
considerably higher than those in the fully developed 
region. This is because the developing boundary layers 
are thinner and offer lower thermal and hydro dynamic 
resistances compared to those for the thick boundary 
layers associated with the fully developed flows. ' Wavy 
and interrupted fins have boundary layers developing after 
each interruption. With a proper design, the resultant 
heat transfer coefficients and heat transfer rates are 
significantly higher at the same pressure drops for wavy 
and interrupted fins compared to those for the plain fins. 
Hence, these fins use the materials very efficiently if 
the design constraints allow the use of such fins. 

In a tube-fin exchanger, tubes of round, rectan- 
gular or elliptical shapes are generally used. Fins 
outside the tubes may be categorized as 

1. Normal fins on individual tubes 

2. Longitudinal fins on individual tubes 

3. Continuous (plain, wavy, inter runted) fins on an 
array of tubes. 

Fins inside the tubes are classified as integral or attached 
f ins . 

3.3 HEAT EXCHANGER DESIGN METHODOLOGY 

Two most common heat exchanger design problems 
are the rating and sizing problems. For an existing heat 
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exchanger, the performance evaluation problem is referred 
to as the rating problem. To arrive at a design of a new 
exchanger to meet the specified performance, is referred 
to as a sizing problem. The objective in rating is either 
to verify Vendor's specifications or to determine the 
performance at off-design conditions. In this problem, 
the following quantities are specif ied: the exchanger 
construction type, flow arrangement, overall core dimen- 
sions, complete details on the material and surface 
geometries on both sides including their heat transfer 
and pressure drop characteristics, fluid flow rates, inlet 
temperatures, and fouling factors. The designer's task 
is to predict the fluid outlet temperatures , heat transfer 
rate, and pressure drop on each side. 

Sizing involves * design of a new heat exchanger 
for specified performance within known constraints and 
minimum objective function. Here the fluid inlet and 
outlet temperatures and flow rates are generally specified 
as well as the pressure drop on each side. The designer 
selects the construction type, flow arrangement, materials 
and surfaces on each side to determine the core dimensions 
to meet the specified heat transfer and pressure drop 
requirements. If the problem is not well posed without 
even minimum necessary info tmat ion, the designer arrives 
at the necessary information based on the discussions 



with the customer, his own experience and engineering 
judgements. 

The package approach suggested by .Shah et al. [ 12 ] 
has been adopted in the present work for the heat exchanger 
design. The complete procedure is given in the flow chart 
form (Fig. 4.1 In this method each possible surface geometry 
and construction type is considered to be an alternative 
design. Thus there may be several independent optimized 
solutions satisfying the problem requirements. Engineering 
judgement, comparison of objective function values, and 
other evaluation criteria are then applied to select a 
final optimum solution for implementation. Thus after 
ariving at the possible alternate designs, the total number 
of constraints for each design are formulated. This in- 
cludes the customer's explicitly specified constraints 
like fixed frontal area, ranges of heat exchanger dimen- 
sions and implicit constraints like minimum heat transf er,Ctn<l 
allowable pressure drop. Starting with a set of dimensions 
for the geometry and input operating conditions, the rating 
problem is solved. The output from the heat exchanger 
calculations is fed to the optimization package where the 
constraints and objective function are evaluated. New 
values for the design variables are subsequently generated 
and heat exchanger calculations are repeated. The itera- 
tions are continued until the objective function is 
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optimized within the desired accuracy and all the cons- 
traints are satisfied. 

A glance at fcteie Fig. 4 indicates two computer 
program packages needed for the optimisation are the heat 
exchanger design programs and the optimization oackage. 

It should, be emphasized that the heat exchanger programs 
must be working properly .before any attempt is made to 
implement the optimization package. Otherwise, " garbage- 
in gospel-.out" will result from the " black box". 



CHAPTER 4 


OPTIMUM DESIGN OF A COMPACT HEAT- EXCHANGER 


4 . 1 INTRODUCTION 

A variety of heat-exchanger design-methods have 
been proposed to determine the best performance of a heat 
exchanger based on specific performance criteria. However, 
little work [ 12 , 14] has been done applying the concepts 
and techniques of non-linear programming to optimizing© 
heat exchanger design. These techniques offer a syste- 
matic approach for improving heat exchanger performance 
by allowing trade-offs to be made between competing objec- 
tives such as pressure drop, power, cost, volume and 
weight. The optimization package developed in Chapter-2 
is quite general and may be used for any design problem. 
The present chapter is devoted to test the applicability 
and efficacy of the general package for the design of a 
compact heat exchanger for a gas-turbine plant. Relevant 
datas for the design parameters have been taken from 
Appendix-B of Kays and London [id]. The objective is to 
minimize pumping power requirements with constraints on 
the exchanger size, volume and the total heat transferred. 
The design variables are the hot and cold fluid capacity 
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flows and linear - dimensions of the heat exchanger. The 
£ — NTU approach is used to evaluate the heat exchanger 
performance. The design procedure outlined in this chapter 
is quite general and suitable for any heat exchanger 
application. 


4.2 NOMENCLATURE 


fra 


f rw 


A 

A 

B 

C 


w 


mm 


max 


pa 


pw 


w 


w 


g. 


Exchanger total heat transfer area on air side 
Free flow area 

Exchanger total fin area on air side 

Exchanger total frontal area on air side 

Exchanger total frontal area on water side 

Exchanger breadth (Fig. 5) 

Flow-stream capacity rate of air 

Flow-stream capacity rate of water 

Minimum of C or C 
a w 

Maximum of C or C 
a w 

Specific heat at constant pressure for air 

Specific heat at constant pressure for water 

Friction factor on air side 

Friction factor on water side 

Air stream mass velocity 

Water stream mass velocity 

Proportionality factor in Newton's second law 
Heat transfer film coefficient on air side 
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h 


w 


H 

k. 

c 

k 

x 

k. 


w 


L 

X f 

M 

m 

P 

a 

P 

w 

Q 

°D 

R 

cn 

r 

a 

r 

w 


h 


U 

V 


v 

w 


Heat transfer film coefficient on water side 
Exchanger height (Fig. 5) 

Thermal conductivity of air 
Thermal conductivity of water 
Thermal conductivity of fin material 
Exchanger fluid-flow length 
Fin length 

Molecular weight of air 
A fin effectiveness parameter 
Pressure on air side 
Pressure on water side 
Total heat transferred 
Total desired heat transfer 
Universal gas constant 
Cana city ratio C . /C 
Hydraulic radius on air side 
Hydraulic radius on water side 
Temperature of cold fluid (water) 

Temperature of hot fluid (air) 

Overall thermal conductance 
Volume of heat exchanger 
Specific volume of air 
Exchanger width (Fig. 5) 
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G reek Symbo 1 s 

a Ratio of total transfer area on one side of the 

exchanger to total volume of the exchanger 
^ Fin thickness 

e Exchanger effectiveness 

^ Fin temperature effectiveness 

>1 Total surface temperature effectiveness 

a Ratio of free flow area to frontal area 

4 Viscosity coefficient 

p Density 


Dimensionless Groupings 

Re Reynolds number (4 r G/ju) 

Pr Prandtl number (u C /k) 

NTU , Number of transfer units (AU/C min ) 
j Colburn factor (h P / G C ) 

ir 

Subscripts 


a Air side 

avg Average 

i Inlet conditions 

o Outlet conditions 
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4.3 PROBLEM FORMULATION 

4.3.1 Design Problem 

The heat exchanger configuration is shown in Fig. 5. 

The inlet pressure, temperature and flow stream capacity 

rate of the hot air are T (126.7°C), P (275.8 kP ) and 

hi a a 

C respectively. The values in the parenthesis have been 

C A <5flS TcR&*rv6 HErffT 

chosen for a specific application^ The hot air is cooled 

by water, pumped at the rate of C and at ambient teraoera- 

w x 

ture of T . (15.6°C) and pressure P, (101.04 kP ) . It is 

required to minimize the total fluid pumping power for this 
exchanger by varying the heat capacity rates C , C on 

3 W 

each side and also varying the height, width and breadth 

of the exchanger. It should be spaced within a volume of 

V and it has a heat transfer rate of atleast CL,, 
max D 

4.3.2 Selection Aspects 

First of all, a cross flow arrangement with both 
fluids unmixed is selected. For a complete optimization, 
every possible combination of available heat transfer sur- 
faces should be considered in order to make, reasonably.- 
sure j that the optimum design is obtained. This may be 
impractical and £<zf® time consuming. For the given pressure, 
temperature ranges and space constraint, as selected for 
the present work in 4.3.1,. a finned-flat tube core confi- 
guration is suitable [ 15 ] . Once the basic surface geometry 
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is selected, the designer gets the values of various 
geometric parameters, as listed helow, that enter the 
heat transfer calculation. 

1. Hydraulic radius, r 

2. Heat transfer area/Total volume, a 

3. Free flow area/Frontal area, a 

4. Fin area/Heat transfer area, A^/A 

5. Fin thickness, $ f 

6. Fin conductivity, k„ 

f 

7. Fin length, 1^ . 

These are to be determined for the hot and cold fluid 
sides of the heat exchanger. The heat core selected in 
this application uses fins only on the air side for the 
obvious reason that the heat transfer film coefficient is 
lower on this side. The friction and colburn factors for 
the flow inside the flattened tubes are obtained from 
[l5] assuming the ratio L/4r ^ = 40, and are given in 
f igures 6 and 7 . 

4.3.3 Optimization Prob lem Definition 

The design variables are exchanger width W, 
height- H, breadth B, and the capacities C a and C^. By 
appropriately scaling and nondimensionalizing these vari- 
ables, a new set of parameters may be defined which have 
the same order of magnitude. This makes the operation 
of optimization procedure more convenient. Defining 



X 


L 


9 


Fig 6, Friction & Colburn factors forAir side 



Fig 7. Friction & Colburn factors for Tube side 
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L 


(V 

max 


) 


1/3 


Z 



( 11 ) 


the variables W, H, B, and are nondimensionalized 
and the design vector is defined as follows: 

X. = c / z 

1 a 


X = C /10Z 

2 w 

(12) 

X, = W/L 


X 4 = B/L 


X„ = H/L 

D 



Representing the pumping power as 0, the minimization 
problem is: 

Minimize 0 ( X ) (13) 


Subject to the following constraints. 


i) 


ii) 


Tc maintain the minimum value of the heat 
transfer rate, Q > 



To space the exchang-er within a -certain volume 


V 

*2 


— ^max 

_ V 

" V 
max 


1 < 0 
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iii) 


Realizing the need for a meaningful configuration, 
the following constraints are imposed: 

Width should be less than twice the breadth, and 

Height should be less than thrice the width 
i -e. 

x 3 < 2 X 4 

X<_ < 3 X 

5—3 


f3_ 

2X ; 

h. 

3X. 


1 < 0 


- 1 < 0 


iv) 


Upper bounds for the capacities are given some 
practical values say BIG 


C < BIG 
a 


C < 10 BIG 
w 


The water capacity rate is taken to be of one 

order higher than air capacity rate. 

X. 


g 5 ~ BIG 


- 1 < 0 


v) 


- 1 < 0 


“ 10 BIG 

Upper bound for the remaining variable breadth 
is given as 1/L so as to satisfy the volume 
constraint always: 


g 7 = l - l < o 
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vi) Non-negativity constraints are finally imposed 

on all the five variables. 

g 3 = - Xj < 0 

g 9 “ “ X 2 < 0 

g 10 = - x 3 < 0 

g ll = “ X 4 < 0 

g l2 = ~ X 5 < 0 

4.4 PERFORMANCE CALCULATION PROCEDURE 
4.4.1 Overvi ew 


The heat exchanger performance is evaluated in 
an iterative manner as shown in Fig. 3. An iterative 
scheme has to be used since the outlet fluid temperatures 
and heat exchanger effectiveness are not known at the 
outset. A value of 0.75 is assumed for the effectiveness 
e at the start of the iteration scheme. The outlet fluid 
temperatures can then be found. The fluid properties vary 
with temperature and are evaluated at the mean of the 
inlet and outlet temperatures. The fluid properties 
appear in the expressions for mass velocities (G), and 
Reynolds (Re) and Prandtl (Pr) numbers. The Colburn ' j 
and friction 'f ' factors are function of Re. With j and 
Pr known , the film coefficient h can be computed, then 
the overall conductance U, and finally the number of 
transfer units NTU. Exchanger effectiveness e is a 



I As sura# Effectiveness 
L £o 0.75 

P ' Calculi Outlet^ 
i r Teraj ergtures _ 

Evaluate Fluid “ 

< , f 

i I roperties ! 


Determine Reynolds 
I & Frandtl Numbers , 

Compute Colburn *r 

; & Friction Factor® j 

L. , J 

! 'Estimate Fils"' 

I 

j Coefficients 

— 1 J 
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function of NTU and capacity ratio ^ 

6 Sre aVallable f °t different flow configurations. The 

computed value of e „ 

e ot e is then compared with the assumed 

value for concurrency. 

4.4.2 Core Dimensions 


The geometrical parameters for the selected 
finned-flat tube surface 11. 32-0.7 37-SR [l5 ] are: 

i ) -Air side: 


r 

a 

= 0.00288 

a 

a 

= 270 

A f /A a 

= 0.845 

a a. 

= 0.78 

6 f 

= O.OOQ33 

k f 

= 100 


= 0.01875 


(14) 


ii ) 


Water side: 


w 


a 


w 


a 


w 


= 3.06 x 10 3 
=42.1 
= 0.129 


(15) 


Prom Fig. 5, the frontal areas for the air and water, sides 
of the intercooler can be readily computed. The frontal 
area is the cross- sectional area normal to the direction 


ms 2 
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of the flow. 


A_ = tajh 
fra 


= X 3 X 5 1/ 


A,. = WB 

f rw 


( 16 ) 


= X 3 X 4 L ‘ 


The total volume of the heat exchanger is : 
V = -WHB 

= X 3 X 4 X 5 L 3 


(17) 


4.4.3 Outlet Fluid Temperatures 

From equation (12) we have 


w 


= x 1 z 


= 10 X 2 Z 


(13) 


Defining C . to he the minimum of C and C , the maximum 
min 3. w 

possible amount of heat that can be transferred by the 
heat exchanger is given by: 


Q 


'max 


= c . (T, . - T .) 

mm hi ci 


(19) 


If the intercooler effectiveness is e * then the actual 
amount of heat exchanged between the hot and cold fluid 
streams is: 


Q = * % 


ax 


e C . (T, . - T . ) 

mm hi ci 


( 20 ) 
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The heat balance equations are; 


Q = C (T, . - T, ) 
a hi ho 


= C (T - T . ) 

W CO Cl 

Substituting (20) in (21), the outlet temperatures 


T = T 
ho hi 


SC. (T, . - T . ) 
mm hi ci 


T = T . + 
co ci 


e C . (T . - T . ) 
mm hi ci 


The optimizer will manipulate the design variables in 
order to find optimum/ ensuring at the same time that the 
above energy balance is satisfied. 


4.4.4 Fluid Properties 


For heat exchanger applications, it is a normal 
practice to estimate the transport properties, which are 
temperature dependent, at the bulb temperature which is 
the average of the fluid inlet and outlet temperatures. 
The bulk temperatures for the air and water streams are: 


cavg 


havg 


T + T . 
co ci 


T 4 - T 

ho + hi 


(23) 


Experimental data for the fluid properties of air and 
water at different temperatures are readily availablef 19] 
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In this optimization study, polynomial expressions are 

-i-itted to the experimental data using the least squares 

technique (Fig. 9). The perfect gas equation, Pv=RT/M, 

is used to determine the specific volume of air at the 

inlet, outlet, and bulk temperatures with the pressures 

assumed to remain constant at P . 

a 


v< = 53.345 (460 + T hj ) 
133 P 

a 


v 

o 


V 

avg 


5 3.345 (460 + T, ') 

hO / 


144 P 


5 3. 345 (460 + T, ) 
havg 

144 P 


(24) 


4.4.5 Reynolds and Prandtl Number 


The mass flow is given by dividing the heat 
capacity C by the specific heat C^. With the free flow 
area known, the mass velocities for air and water are: 


G 

a 


C A 
pa ca 


C a 
pa u a fra 


w 


w 


C A 
pw cw 


C CT 
pw w f rw 


(25) 


The Reynolds and Prandtl numbers can now be computed as: 



Re 


Re 


w 


Pr 


Pr 


w 





JLt C 
w pw 


w 




(26) 


4.4.6 Fricti on and Colburn Factors 

Polynomial expressions. Fig. 9 are fitted to the 
experimental data as shown in Figures 6 and 7 by least 
square method. In practice, the experimental data is 
approximated by several polynomials which are valid for 
a small range of Reynolds number. Since most compact heat 
exchangers operate in the transition flow regime (100 < Re 
< 10,000) it is important to get good polynomial fits for 
f and j in this flow regime. 


4.4.7 Heat Transfer Film Coefficient 


Once the colburn factors for the hot and cold 
fluids are known, the film coefficients are calculated 
from the following equations; 

, _ J ’a G a C pa 

a ~ (Pr ) 2/3 
a 

j G C 
J w w pw 

“ " (Pr ) 2/ 3 
w 


( 27 ) 
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4.4*8 Fin Ef f ectiveness 


For finned-tube compact surfaces, the fins are 
assumed to be straight uniform fins. The fin effective- 
ness for this type of fin is shown to be [ 20 ] . 


Hf 


where 


m 


tanh (ml f ) 
~mi~ ~ 


2h 

= — S' ) 

f f 


1/2 


(23) 


The surface effectiveness of the heat transfer surface 
on the air side is given by [ 20 ] ; 


n = i - ~ (i -t) ) 

ri r 

a 


4.4.9 Exchanger Effectiveness 


(29) 


The conductance of the heat transfer surface can 
be computed from the film coefficients and the surface 
effectiveness. In computing the overall conductance, the 
wall resistance of the surface is neglected. 


1_ 

U 

a 



a 

_v 

a 


w 


w 


( 30) 


a 

Equation (30) gives the overall conductance based on the 
air side heat transfer area. The number of transfer 
units (NTU) is 


NTU 


A U 
a a 


mm 


a V U 
a a 


C 

mm 


(31) 
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Let C max re P res ent the maximum C and C . Then the 

a w 

capacity ratio c* is; 

„ C . 
bt- mxn 

C _ - ( 32 ) 
max 

With NTU and C* known, the effectiveness e for the cross 
flow intercooler with both fluids unmixed, is determined 
from the series solution of Mason [ 2 1 ] 


1 £ [ l-e~ NTU 2 (MTU ) m 

(NTU ) C* n=o m=o m 


- -(NTU)C* J 
l-e L 

m=o 



(33) 


This series is rapidly convergent. The value of e from 
(33) is compared with the value assumed at the beginning 
of the iteration. If these two values differ by more 
than 0.01, the performance calculations are repeated, 
using the new value e for the next iteration. A few 
iterations are needed to obtain agreement between the 
assumed and computed effectiveness of the compact heat 
exchanger. 


4.4.10 Pressure Drops 

The pressure loss in the fluid streams is con- 
tributed by t)?®. change in fluid momentum, resulting from 
the changes in the cross-sectional area at inlet and 
exit, and by viscous friction. The pressure loss for 
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the air stream is given by [ 15 ] 

Ap 


= [( 1 + a a b £ -l) + f a ^^a] 

1 Cl 


c a 


From the definition of hydraulic diameter: 


A c ' r a “ r a 


(34) 


(35) 


On substituting (35) in (34), we get 


AP 


G 

< 

3600 


(- .3 ) 2 v 

v n/ron ' v * 


a 2g P 
^c a 


_ 0 v f XL v 

[(1+ ® a 2 ) c^-n ] 

i a i 

( 36 > 

The pressure loss associated with the water stream is 
mainly due to surface friction and is given by: 


ap = (-§—-)' 

w 3600 


f H 
w 


2 g p r 
w 


( 37 ) 


G _ f X C L 
, w n 2 w 5 

“ ^3600 ‘ 2 g p r 


'c w 

where, g c is the gravitational constant, and p is the 
average mass density of water. 

The pumping power for each media is the product 
of pressure loss and volume flow rate, and the total 
power is given by: 


C . V 

_ a avg 

0 = AP a ’ C 

pa 


ap c 

w w 

P C 
pw 


(38 ) 
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4.5 SOLUTION TECHNIQUES 

Three dirferent constrained minimization tech— 
niques were applied to minimize the pumping power (Eqn. 
13). They are Interior penalty function. Exterior pena- 
lty function and the Method of feasible directions. For 
penalty function methods, the multidimensional minimiza- 
tion schemes tried were Steepest descent, Powell's conju- 
gate directions, Fletcher and Reeves' conjugate gradient 
and the Davidon-Fletcher-Powell variable metric method. 

For obvious reasons. Univariate method was ignored. These 
methods were made to select the Quadratic interpolation 
procedure for the unidimensional minimization. A graphics 
program, using the GPGS (General Purpose Graphics System) 
routines, displays on the graphics terminal the heat 
exchanger configuration according to the optimized dimen- 
sions. This may be used '.also to give the initial solu- 
tion for linear dimensions of the exchanger, at the 
start of the execution of the optimizer, to suit the 
designer's taste for configuration. 

The heat exchanger performance analysis is put 
in g separate routine. Appendix— A. Finite difference 
schemes were utilized to get the gradients for various 
search techniques in the optimizer. All the input para- 
meters used for the heat exchanger performance analysis 
have been defined in sections 4.2 and 4.3. 
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A sample of the input data, for the optimizer, used to 
define the variables and initialize the flags is given 
here* Results are given in the graphical form. Varia- 
tion of pumping power and eff ectiveness with the volume, 
has been plotted in figures 10 and 11, and a detailed 
discussion is reserved for the next chapter. 


List 

of Input Data 

Significance 

NN 

= 

5 

Number of design variables 

NC 

S 

12 

Total number of constraints 

NI 

= 

12 

Number of inequality constraints 

X( 1) 

= 

2.5 

Starting point (Feasibility 

X ( 2 ) 

= 

1.0 

is a must for Interior penalty 

X( 3) 

5= 

0.5 

and Feasible direction 

X( 4 ) 

= 

0.3 

methods ) 

X( 5 ) 

rr 

1.25 


VMAX 

rr 

16.0 

Max. volume allowed in a specific 
run 

BIG 

rr 

5.0 

Max. capacity rate of both fluids 

TT 

= 

0.1 

Initial, step size 

TMAX 

rr 

0.5 

Max. step size allowed 

LIM 

rr 

15 

Max. number of iterations permitted 

FF 

= 

0.01 

Percentage change in variables 
for FEM computation of gradients 

ACC 

= 

1.E-04 

Accuracy desired 

R 

rr 

10 

Initial, penalty factor 
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FC = 0.1 

= 10 
W(1:NC)= 0 

NPERT = 1 


Factor decreasing R for IPENAL 

Factor increasing R for EXPEN 

Push off factors set to 'O' for 
linear constraints 

Ships the perturbation scheme 



CHAPTER 


5 


RESULTS AND DISCUSSION 

5 . 1 TRADE-OFF STUDIES 

The compact heat exchanger optimized in the 
last chapter has the same configuration as considered 
by Kays and Londan in example-1 [l5]. The objective of 
the study is to find what reductions in pumping power 
or volume is possible for this design, while still keep- 
ing the same heat transfer performance. The heat exchan- 
ger is designed by varying it's dimensions and the mass 
flow rates of the two fluids viz, air and water. The 
results, as given in figures 10 and 11, are indicative 
of the characteristics of a single optimization. Fig. 

10 indicates the variation of pumping power with the 
exchanger volume and Fig. 11 is drawn for effectiveness 
of the exchanger versus volume. Each point on the curve 
in Fig. 10 corresponds to the optimum! value of the pump- 
ing power for a specific size of the exchanger. Simi- 
larly, in Fig. 11, the optimum performance of the exchanger 
is indicated for a particular volume. 

However, an expanded optimization study can 
easily be done by changing the appropriate constraints. 
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For example, if one is interested in finding the iluid 
pumping power requirement as a function of the exchanger 
total volume, v max is varied and the resultant optimize 
tion problem is solved again by the methods used in the 
last chapter. Optimization results of eqns.(l3) for 
different values of showed that the constraints on 

linear dimensions of the exchanger are always satisfied 
and hence, x, , x 4 and x 5 are determined by solving those 

constraint equation as: 

x^ = 0.8736 
= 0.4368 
x 5 = 2.6207 

. ( 13 ) is now reduced to a 2 

The minimization problem xn 

variables problem and can be restated as 

0 (x^ ' x 2 ^ 


Minimize 

X 

subject to 


g l (x l ' V " 1 


Q (x ^ • X 2 ^ 


, .* - solved by the same three 

This two-variables problem is sox 

.. refrained minimization for diffe- 
different methods of con^-tr- 

•hr-ade— of f curve between pumping 

rent values of A 

is ©resented in Fig. 12, 
power and exchanger volume - 

_ w ith volume in Fig. 13, and 

variation of effectiveness *1 - 

fln „ rpte s with the volume is shown 

variation of mass flow r 

in Fig. 14. 
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It Is obvious from Fig. 12 that the pumping power 
increases rapidly and the effectiveness decreases slowly 
as the volume decreases. The asterisk in Fig. 12 repre- 
sents the pumping power and the volume corresponding to 
the unoptimized exchanger studied in reference [l5j. For 
the same volume/ the optimum design reduces the power 0 
from 31.02 HP to 9.03 HP. If however, the pumping power 
requirement is kept the same, say 31.02 HP, the volume of 
the optimized exchanger is only 0.1614 rrT as compered to 
the 0.4785 m“, for the unoptimized case. Fig. 14, indi- 
cates the behaviour of the flow rates with volume of the 
heat exchanger. It is interesting to note that there is 
only a slight variation in the magnitude of the flow rates 
of both the fluids with volume. 

5.2 A COMPARISON OF OPTIMIZATION METHODS 

With the large number of optimization strategies 
available, the natural question arises - " which Is the 
best strategy?" -At present, a comparitive study of opti- 
mization methods applied to heat exchanger problems, in 
particular, is not available. However, several studies 
have been performed to evaluate the effectiveness of 
various optimization techniques by applying them to spe- 
cially prepared test problems given by Himmelblau [ 2 ]. 
Most heat exchanger optimization problems are constrained 
the following observations are presented for 


and, hence. 
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the constrained methods only. 

(i) The sequential penalty approach and a gradient 
based search are more efficient when a large- 
number of design variables are involved. The 
sequential exterior and interior penalty function 
methods work well with the gradient based search 
procedures even when using finite difference 
gradients. Heat exchanger performance is not 
generally described by a differentiable function 
and, therefore, internal finite difference schemes 
are necessary. The exterior penalty function 
method does not require a starting point which 
satisfies all the constraints. This is a distinct 
advantage in starting a heat exchanger design 
optimization for an unfamiliar application. 

(ii) The method of feasible directions is rapid and 
effective in solving smaller problems, with 
perhaps up-to three or four variables and having 
mostly linear constraints. This method, however, 
may have significant difficulty when a large 
number of desiqn variables or nonlinear -cons - * 
traints are involved. Constraints on pressure 
drop heat transfer rate or weight, for example, 
are nonlinear. Geometric constraints are gene- 
rally linear. 
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5.3 CONCLUSIONS 

1. Well structured optimization program.s-package 
has been developed and tested for it's reliability, 
flexibility and accuracy with standard functions, avai- 
lable in the onen literature [ 1,2,3], 

2 . A direct transfer compact heat exchanger for use 
in gas turbines is designed for minimum pumping power 
subject to constraints on volume, dimensions, and heat 
transferred. The minimization problem, formulated in 
(13) in Chapter 4, is guite general and can be used for 
any heat transfer application. The e-NTU approach is 
used for evaluating the heat exchanger performance. 
Polynomial expressions are found to approximate the 
experimental data for the fluid transport properties and 
the friction and colbum factors for the heat transfer 
surfaces . 

3. The optimization study of the compact heat 
exchanger provides a design which teas got much superior 
performance. 
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APPENDIX - A 


i USER'S MANUAL 




LIST OF SUBROUTINES 


SUBROUTINES 


PURPOSE 


UNCONSTRAINED' OPTIMIZATION : = 


UN IV 

STEEP 

CQNDIR 

CQNGRA 

DFPH 

QUAD 

GOLD 


Univariate search 

Stefepest descent search 

Can jugate directions search 

t Powell's Method 3 

Conjugate (gradient search 

T Fletcher - Reeves Method 3 

Variable Metric Method 

C Davidon * Fletcher - Powell MeltiJod* 3 

i -■ D minimisation oy Quadratic 

interpolatoin scheme 

1 D minimisation py Golden section 

method 


UNSTRAINED OPTIMIZATION ;* 


IPENAL 

EXPEN 

FEAS 

START 

LINEAR PROGRAMMING. S = 
linear 

INTERNAL SUBROUTINES := 


CONSTR 
GRAD » 

MODIFY 

COMVRG 


IRAQI 


Interior penalty functions net odd 1 
Exterior penalty functions net ndd- 
Me.tnod of f easiole directions 
Gets a feasible starting ooint 


Linear programming oy simplex? neitnbd- 


Deternines constrains violated 

Pacfcage routines for determing gradient by 

forward differences scheme: 

Modifies the objective combining thtei 
constraints „ , i 

Cnetics for convergence oy perturbation schemf 


USER WRITTEN SUBROUTINES J= 


FUNCT 

GRADU 


User written function subroutine] 
User written gradient 1 subroutine! 



COMMON STATEMENTS 


- Sole common statements have to oe incl ude’d- in user 
written orogramies 

1. ' C0MMON/A£jWAYS/TMAJC,METHOD,KFEAS,LIM/ NC,NI , I3DT, MPERT 

Mast oe Included in the mainline. 

2. commqn/penfun/r,uf 

Insert in mainline program and function subroutine 
F’JNCI wh'sn penalty functions are used 


3,- 


CDMMOM/CQN VEC/G 
DIMENSION G(?5) 


These statements must be included in function s 
sabrouitne FUNCT if constraints are involved-.' Tn'ei diaensi 
of G. is al*a vs set to 25 in order to overcoiei tnii problea 
of transferring the value of G from FUNCT to internal 
subroutines. 


4.* COMMON /COM AT/DG 
DIMENSION DG(-,-D. 

Placed in mainline program and used with: gradient . 
subroutine GRADU if constraint gradients are’ nete’ded as in 
subroutines START, FEAS .. 

, . Matrix DG must be properly dimensioned using a 
appropriate integers, it's a matrix .of order CMC', NN) , tfi 
elements of *nica represent the components of tn’ei 
constraints gradient. 



EX ’URGES OF STANDARD TEST FUNCTIONS 


TEST FUNCTION IS TAKEN FROM FOX r Ref 1 3 Page 89 
COMMON/ ADNAXS/TMAX, METHOD, KFEAS,LIM,NC, HI f T3UI f XPERT 

com'Um/men/ff.accJgrad 

DIMENSION X C 25 

OPEN (UNIT*! , FILE.= * NN. " ) 

NN = 2 
NC = 0 
U = 5 
X{1)=0. 

X(2)=Q. 

FORMAT! *X( I) := %9£IPE9.2)) 

41 If I. Ms 2 
TMAX=.5 
n=.i 

ACCM.G-05 

HPEftFs! 

LIM=15 

METHOD®! 

Ffs.n, 

F >r * OHi 'Sm» f 

TO USE SRAOU , IGR5D = 2 

room 

XFSASs! 


GAOL- UH1V(4INIM.X,TT,NN,F) 
rfRirpt ,io)(xci)a=i,2).F 
FORMAT! / r X 1 := J ,2UPE 


CLOSE! UNI T=l) 

STOP 

END 


2CIPE14.6)//' t F 1 


1PE14.6) 


SUBROUTINE functc X,F) 

DIMENSION X(2) 

F=t5.*X(D-4-20.*XC2m(l)*2 + 10.*XC2) A 2 + X(i)*2-2',*X(m5. 

return 

end 

SUBROUTINE GRAD J!X,S, SUM) 

OIMENSTON X(2).S!2) 

S(U=-4D.*X(1 )*3 + 40*Xl2)*X(n-2.*X(l)+2. 
S(?J=2D„*X(l) A 2-23.+X(?) 

5U*=Stt) A 2+S(2) A 2 
TFCSJvjOl ,2,1 
SJ« = 1 . 

RETURN 

END 

RETURN 

END 




=. FOX l Re* S 1 3 PROBLEM F5 PP 89 


SOURCE **• L Ket 11 J 

HO OF VARIABLES »=' TWO 

CONSTRAINTS «*■ 

TEST FUNCTION S = 

minimize: 

P(X) = io.«$tl)-4 ;T??5*J t 5 , o X<1) ‘ 2 + 10 -* xt21 2 

v * X( 1 ) 2 **. 2.*X(1) + 5,0 

. tutc smim-tION belongs to rosembrock's banana 

NOTE * THIS F VALLEY FUNCTION FAMILY 

STARTING. POINT S* 

t % r« :=• C 0«0 ' °*° 1 
RESULT GIVEN S s 

' J " {=: Siii ’ *• 1 

RESULTS b3TAOEOC = 

3Y UNIV | , s r 9 3999 E-01 # 1.0000 E + 90 1 

f l 1 l 4.5ooo £+00 

BY STEEP. r 9.3890 E-01 , 9.7736 E-01 J 

C } 1 \l 4.0001 E+00 

BY CONOIR r 9.3990 E-01 , 9.9981 E-01 1 

1 J 1 |s 4.0000 E+00 


l: 1 <U ' *• 1 


9 

V 


BY CONGRA ^ J 
F 

BY DFPM . „ , 



3.1308 

.0355 


E-01 

E + 00 



45S2 E-01 
030 E+00 


b.534l E-01 J 
8.9172 E-01 1 


so(JRCE 1 = D.M.HIMMELBLAU C Ref i 2 1 PROBLEM S3 25 PP 427 

SO OF VARIABLES : = FOUR 

CONSTRAINTS' ?=• NOSE 

test function ;«• 


minimize: 

F(X) 




*2 * 
*4 


STARTIN 3 . POINT 
t X 3 


| RESULT 


STYE'S 


c X J« 
F 


:= £- 3.3 
: * 

; a. [ D ,3 
:=. 0.0 


- 1.0 , 0.0 , 1.0 1 
0.0 , 0.0 , 0.0 1 


RESULTS OBTAINED** 
BY STEEP- .. 

f * a 

BY CONOIR r „ 

f X |i s* 

F S * 

BY CDNGRA* 

[XI< 

■’ If 1 * -* 

BY DFPM „ 


r 1.544E-01,-1.617E-O2,7.949E-O2,8.235£H>2 1 

1.I44S-03 

[ 3 , 4265 - 02 ,- 3 . 472 E- 03 f 3 . 635 E- 02 , 3 . 7 39 SKJ 2 1 
3 . 4322-05 

[ 1 , 0482-02 #-l .038E-Q? .5.859E-Q2 #5 . 943E, -02 1 

3.114E-04 

r i.535E-02,-1.661E-03,-i.457E-02,-l .413E-02 3 
1 .056E-05 


THIS PROBSLM IS TAKEN FROM S. M.HIMMEL8LAU CRef i2J N3:24 PP42 


IGRAQU 


CDMM3N/ALrfAYS/TMAX,MErHOD,KFKAS,HM,MC,Ml,I3Ur,NP£.Rr 

COMMOM/PEMFUM/R^F 

:om^3«/mew/ff,a:c, igrad 

CD M'M 3 N / CD M A I/DG 

DIMENSIDM X C 2 3 ,DG(2, 2) 

OPEN (UNIT*!, FILE* '»(!.*) 

MN = 2 
NC = 2 
MI = 2 

XCt )=2.0 
X ( 2 ) =2 , 0 

rfRlTEC! .5)CX(I)#I=1,MM) 

FORMAT! *X(I3 := , 9 ( IPE9 , 2) ) 

MI MIMs 2 
FMAXs.5 

Yf~ 1 

ACC = 1.5-!>5 
MPS'RTtei 
L IM= lb 
METii33=5 
FF=. Dl 
T‘GRA.3®! 

TD iiSE 3RADU , I GRAD If = 2 

TD«r*S 

KFSA5=2 

FC=.i 

CALL- IPEMALI MINIM, X,TT,NN,FC,F) 

WRtTS’C 1 # 1 3) ( XC I } , I s ! , NN) , F % . - 

FORMAT!//' T X ] := * , 3 ( l PEI 4.6) // * ( F 1 := ',1?S14.6) 
CL3SS'umr=i) 

STOP 

END 

SUB ROUT I NS FfJNC T ( X , F) 

CDMMOM/PEMFUM/R,UF 
CDMMOM/CDMVeC/G 
OIMSMStOM X(3) ,3(25) 

»-=(*( t)-2.)~2 + CXC2)-t .)*? 

3U)=X(1)*2 - X ( 2 ) 

G(2)=X(l) + X(2 ) - 2. 

’TF»F 

return 

END 

SUSR3UTIN2 GRAD J( X, S, SUM) 

CDMM3N/CDMAT/DG 
DIMENSION XC2) ,S(2) ,OG(2,2) 

S(l)=-2.*(X(l)-2.) 

S(2)=-2.*(XC2)«1.) 

SUM=S(t)*2+SC2)*2 
IF(SJM)1 ,2,1 
SUM=t . 

33(1 » 1 ) = 1 . * X( 1 ) 

3G(1 ,2)=-l . 

3 u ( 2 , 1 )*1 . 

33(2 ,2) = 1 . 

return 

EM 3 


SOURCE 


JS O.M. HIMMELBLAU [ Ref ! 2 } PROBLEM S3 2 % PP 426 


HO OF VARIABLES : = TWO 
CONSTRAINTS := 


LI SEAR INEQUALITY S ONE 

N0N LI HEAR INEQUALITY S ONE 


TEST FUNCTION 
MINIMIZE 

PCX) 

311 ) 

3.(2) 

STARTING POINT 
[ X 1 
RESULT GIVEN 


* (X( 1 ) -2. ) *2 + (X(2)-l.)"2 

* X(l)*2 - XC2) <3 0.0 

= XU) t XC2) - 2. <* 0.0 


:*■ l 2.0 , 2.0 ] 


CXI :=• [ l. * i. I 


RESULTS OBTAINED:* 

DFPM IS USED' IN ALL THE CASES 
BY IPENAL 1 

r X 1 := [ 1.0001 E + 00 , 1.0005 E+00 1 

F := 9.9958 £-01 

BY EXPEN 

f X 1< := [ 1.0001 5+00 , 1.0005 E + 00 3 

F := 9 .9963 E-01 

T X 1 :* C 9,9795 5-01 . 9.9795 5-31 3 

F := 1 .0041 5+00 


BY FEAS 


i=. s.S.RAO t Ref : 3 3 PROBLEM 7.7 PP 391 




SOURCE 

flO OF VARIABLES := THREE 

CONSTRAINTS **■ 

linear inequality : one 

NON LINEAR INEQUALITY S TWO 
BOON OS ON INDEPENDENT ? THpg;E 

variables 


TEST FUNCTION 
MINIMIZE 

F(X) 

NOTE 


L 3 7 

si 4) 
"f53 

j'l J { 

3<S) 


- x{ 1 ) * 3 - 6.*X(1)“2 + It .*X(1) + XC3) 

: XC2) IS NOT INCLUDED IN FIX) OEFENITION 

: ... 

^ SjiAH.r °-° 

3 -x( 2) <= o.o 

* -X(3) < = 0.0 


STARTIN3. POINT ; = 

C X 1' := t 0.1 r 0.1 * 3.0 3 

RESULT 31 YEN S 3 

r X 1 : s ( 3.24BE-08 , 1.41421 , 1.41421 3 

1 f 1 ;* 1 1 . 4l 422E + 00 

RESULTS OBTAINED: 3 

DFPM is USE O' IN ALL THE CASES 

BY IPtNAL - j .30E-05 » 1 .41465E+00 , 1 ,4l467Ef 09 J 

r • x 1»447B21E+vO 

BY EXPEN 5.52&9E-04 , 1 ,401 93E-02 , 1 .99367E+3D 3 

F := 1.993S8E+00 


BY FEAS 

1 5 1 


;= £ 0 , OOOOE+OO . 1.41422E+00 , 1 .414222*30 3 
:= 1.I1422E+00 


SOURCE I* D.M.HIMMELBLAU t Ref i 2 1 PROBLEM NO 11 PP 406 

NO OF VARIABLES :=■ FIVE 
CONSTRAINTS i= 

LINEAR INEQUALITY : 6 

SOU NOS ON INDEPENDENT 

: 10 

VARIABLES 

TEST FUNCTION' 

MINIMIZE 

PCX) * 5.3573547*X(3)*2 + ,8356891*X(1 5*XC5) ♦37.293239*X 

- 45792.141 

NOTE S XC 23 & XC 43 AREN'T INCLUDED' IN FIX) D£F SNITION 

CONSTARI NTS ARE' ! = 

n.O €* 35.334407 + 0 ,0056858*XC2) *X( 5) ♦ 0.00052S2*XC1)*X£43 
0lO022O53*X(3)*X(53 <= 92.0 

90.0 <* 80.51249 + 0.0071317*Xf 2)*X(5) t 0.0029955*X( 1 ) *XC2) 

0.0021813*X(3)"2 <s 110.0 

20.0 <a 9.300961 +• 0,0047026*X(3) *X(5) + 0.0012547*X(1)*X(3) 

0.00190B5*X($)*XC4) <= 25.0 


73.0 

33.0 

27.0 
27.0 

27.0 


X(l) 

xti? 

X(4) 

XC5) 


102.0 

45.0 

45.0 

45.0 

45.0 


STARTING POINT 
r X 1 


*. [ 78.62 , 33.44 , 31.07 , 44.18 , 3 5.22 ] 


RESULT GIVEN 

C X I 
F 


= [ 78.003,29.995,45.000,36.776 1 
= - 3 0665 . 5 


RESULTS OBTAINED 


BY DFPM 

t X 1 
F 


;= £ 78.595,33.317,30.739,44.162,35.129 1 

;= - 35493.3 


A. 4 COMPACT heat exchanger problem 


COMMON/ALWA YS/TMAX, METHOD, *FEAS,LIM,WC, NT, TOUT, NPERT 

COMMON/PENFUN/R , UF 

COMMON/NEW/FF.ACC, IGRAD 

COMMON /C ALL /NC ALL , VMftX , BIG, VOL 

COMMON/PARA/YY 

COMMON /EXCHGR/EF 

DIMENSION YYC5) 

OPEN (UNITsl .FILES "R. # 3 
CALL RTIWECHI) 

NCftLLaO 

NNr5 

NC*12 

NIsl2 

YYri)r5.5 

YYf2)=1. 

YY(3 )*,§ 

YY(4)». 3 
YY(5)=1.25 
VMAX«14.0 
BIG=5. 


WRITE Cl *5)fYY(I),I«i,53 „ ^ % 
FORMAxC ' t X 1 *s ,9( 1PE9. 2 ) 3 


WRITE (I #6) RIO, VMftX 

FORMAT( f MAX f6lW RATE s« ',F5.2, ' V max I* ',F5.?) 

MINIMs2 
TMAXs.5 
TT*. t 

ACC=t.E-05 

NPERTsI 

LIM»t5 

METHOD=5 

FFs.Ol 

IGRADsl 

TO!!T=6 

KFEAS=1 

PC=iO. 

OsiO. 

CALL EXPEN (MINIM. YY,TT,NN,FC,F3 
WRTTEC 1 , 11) EF, NCALL 

FORMATC 'EFFECTIVENESS ' , 1PE14,6/'N0 OF FUNCTION CALLS ',T53 
WRTTE (1,12 3 VOL 

FORM AT ( 'VOLUME := ',1PF14.6) 

CALL RTIME(MP) 

MTxME-MI 

WRTTFCl .103MT _ „ 

formrtc 'CPU TIME t* ',15, 'ra.sec'/) 

CLOSE(UNIT=l) 

TYPE 100 

FORMATC' WANT GRAPHIC option KY/N3 '$) 

ACCEPT 101, ANS 
IF(ANS.HE.*Y')GOTO 102 
CALL GRAPH 

YOU HAVE TO ASSIGN T T Y 5 FOR GRAPHIC OPTION 

STOP 

END 


i= ',15, 'ra.sec'/) 



c 

c 

c 


c 

c 


SUBROUTINE FIJNCTCYY,F) 

HEAT EXCHANGER PERFORMACE ANALYSTS ROTffINE DEFINING THE 
OBJECTIVE FUNCTION 
COMMON /PENFUN/RR , UE 
COMMON/CONVEC/G 

cqmmon/call/ncall,vmax,big,vol 

COMMON/EXCHGR/EF 
DIMENSION YYC5),G(25) 

REAL NUC,NTU,L 
NCALL=NCALL+t 


DO 5 1*1,5 
IFfYY(I)5lO0,100,5 


ONTINUE 
T HI*260. 

TCTs^O. 

PA*39.1 
P*®14.7 
OD=8,E+06 
SURFACE CHAR 
RH*. 00288 
ALFAH*270, 

SIGMAH*,78 
RC®. 00306 
ALEAC®*2,1 
SIGKAC*,|29 
FIN CHAR 
FA*, 845 
DEL®, 00033 
FL*. 01875 
FKslOO. 

NON DIMEN 
L®VMAX*(l./3.) 

ZsOD/tTHT-TCT) 

SOLUTION 

AF®H®YY(3)*YY(5)*L a 2, 

AFpC=YYC3)*YY£A)*L a 2. 

VOL=YYC3)*YY(4)*YYC5)*VMAX 
El®, 75 
CH*YY(n*Z 
CC®YYC2)*Z*10. 

CMTNsAMlNlfCH.CC) 

CMAX=AMAX1£CH,CC) 

E®E1 

Q=CMIN*E*(THI-TCI) 

TH0®THT-E*0/CH 

TC0®TCT+F*0/CC 

thavg*cthi+tho)/2. 

TC& VCs f TCT + Tmi /*> 

CALL PROPC(TCAVG,'$HDC,CONDC,CPC,nc) 

CALL PROPHfTHAVG,RHOH,CONDH.CpH,UH) 

Vl=53.345*( 460,+THI) / (144, *PA) 
V2=53.345*(460.+THO)/(144,*PA) 
VM=53.345*(460.+THAVG)/(i44.*PA) 

GH=CH/fCPH*SIGMAH»AFRH) 

GC=CC/CCPC*STGMAC*AFRC3 

REC=GC*4.*RC/UC 

REH=GH*4.*RH/UH 

PR^=HC»CPC/C0LDC 

PRH®UH*CPH/CONDH 

DETERMINE COLBURN an* FRICTION FACTOR ERO.M POLYNOMIAL FITS 
Y=ALOGiOCREH*.DOt) 

FH*«. 11 9O8*Y-4.-.044fi28*Y*3.+.26677*Y-?.-.45865*Y+. 58734 
C0H=-.l I873*Y“4. +. 03586 1*Y*3,+.2?707*Y-2,-,478R4*Y-. 0062732 
FH=10.'*CFH)*,01 

coh=iq »ccoh5*,oi 

HH®COH*GH ACPH/PRH* ( 2, /3 . ) 

Y®ALOG10(REC*.0015 
IFCY.LE.O. >GOTO 1 
IF(Y,GE,l,3GOin 2 

FC®-i.2§3§*Y'‘4. + 2.1 64 3*Y-3.-.56969*Y“2.-.87268*Y+. 38628 

CQC=-1.343t*Y-4.+l.603?*Y*3,+.55494*Y-?.-1.2064*Y«. 083788 

GOTO 3 

EC=.3802-.8745*Y 

C0C®-,09691-Y 
GOTO 3 


** * -jy m 



100 


CDC=10 *ccnc)*.oi 

HCsCOC*GC*CPC/PRC*C2./3.) 

FMsSORT£2.*HH/fFK*DEL)5 

YsFM*FL 

FI*iEFF*TANH(Y)/Y 
FFFssf f i »FIMEFF5 

UCONOsEFF*HH$(ALFAC/ALFAH)*HC/CALFAC/ALFAH*HC+EFF*HH) 

NIU=ALFAH»V0L*UC0ND/CMIN 

RsCHIN/CMAX 

CALL EFFECT(F1,NTU,R) 

IF( APS (E-El ).GT..0i) GOTO 9 

PHr?#2!5vi*CU,+SIGMAH*2,)^£V2/Yl-i.)+FH*YY(4)*L»CVM/Vl)/RH) 
1 /64,4 

P§ilRi?PH?CH*i«/CPH+PC*CC/(RHO^IcpC))/( 5 * 50 .*S 6 O 0 .) 

FsPOWER 

UFsF 

EF*E1 

GC1 )*-YY(l) 

G(2>=YY(1)/BTG-1. 

1(4)»YY£2)/C10.*BIG)-1. 

!(5)*-YY(3) 

|=-YYf41 
>*-YY(5) 

,.))-YY{5)/(3.*YY(3))"l» 

11)=V0L/VMAX-1, 

M 2 )=i.-a/G 0 

iTUPK 
ENO 




SUBROUTINE PR0PC(T,R,C,CP,U2 
c POLYNOMIAL fits for transport properties of air 

XsT/100. 

CPs, 001 7822*X*4.-,0l7893*X - 3.+.O67628*X*2-.091686*X+l. 0348 

R*l6,*C-.Q019308*X'*4. + .02l8*X-3.-.13*X*2. + ,080173*X+6.2283) 

U=.045473*X'‘4.- ? 56414*X-3.+2.6)07*X-2.-5.5424*X+5.1444 

C=, 1 *(-.O0332l9*X < ‘4. + .O437O6*X'‘3.-, 29151*X*2+. 88891 *X+ 3. 0051) 

RETURN 

END 

SUBROUTINE PRQPHCf ,R,C,CP,U} 

C POLYNOIMIAL FITS FOR TRANSPORT PROPERTIES OF WATER 

X=T/100. 

Us5.1*C“» 0057 09 13*X*4* + . 05755772 * X* 3. -,?2800553*X*2 + 
t +,6754148*X+. 46343288) 

CP«.l»(.03306576*X-4.-.l9582903*X-3+.3068203R*X*2. + 

1 + ,20886326*X + 2.36525726) 

Cs.01*C.OOf34717*X-4.-.04445982*X-3.-.10930872*X A 2. + 

I ♦ 2.0579495*X +1.48826838) 

RETURN 

END 

SUBROUTINE EFFECT(E,NTU,R) 

C EFFECTIVENESS ^OR CROSS-FLOW EXCHANGER WITH BOTH FLUIDS UNFIXED 

REAL NTU 

SUM® ( 1 . -EXP C -NTU ) ) * ( 1 . -EXP C -«TU*R ) ) 

SUMOsSUM 

N=1 

3 SU«1=1. 

SU«2=1. 

SM=1. 

DO 1 J=1,N 
SM=S«*.J 

SUMl=SnMl+CNTU)*J/SM 
SU M 2=SUM2+(NTU*R) * J/SK 
t CONTINUE 

SUM=vSUM + Cl.-EXP(-NTU)*SUMl)*Cl.-EXPC-NTU*R)*SUW2) 

IF f APS (SUM-SUMO). LT.l.E-05) GOTO ? 

N=»J + 1 
SUMOsSUM 
GOTO 3 

2 E»SUM/(NTU*R) 

RETURN 

END 



c 

c 

c 


5 

6 


7 

40 

15 


n 


32 


12 


SUBROUTINE graph 


THTS ROUTINE USES GPGS * 
COMMON /PARA/YY 
COMMON/TRANS/W1,W2,W3 
COMMON /EXCGHR/EFF 


outi" eS 


for GR1PHIC RUR^ SES 


wl 



. W2C6) ,W3f6) # 


yy(5) 


DATA 

DATA 

DATA 

DATA 


/ *1 o . » iooo 

¥ J/l • 1 1<3 1 1 • > 1 • » 1« ' < * 1 o f 900, / 

W2/-10.,900.,-10. ,9^$0 w" 1 

tf 3/901 ,,1204,, -10, » * u 
cH=yyci)*40O00. 
cc*yy( 2 )* 400 ooo. 

XsYY f 3) 

Y-YYC5) 

ZaYYf4) 

CALL NTTDEYCS) 

CLRDEVf 5) 

VP0RT3 i VI ) 

WINDW3CW1) /N) '§) 

FQRMATC6 ( /) ' W A NT taPTT 0 
ACCEPT 6 ,CAP 
FORMATCA1 ) 

TF(CAP.NF.'Y')COTO 7 
CALL Cf.RDEV ( s } 

CALI# TITLE rt . ?T TNUING *! 

TYPE 40 r . nP CO MT ‘ 

FORMAT ( 5 f / 1 ' ! I RETURN CO K 
ACCEPT IS ,^TV 
FQRMATf A4) m n0 

IFCTTY.EO. 'STOP') GOTO 1 0 


./ 


CALL 

CALL 

CALL 

TYPE 


'$1 


CLRDEV ( 5 ) 
MARGIN 
VPQRT3 (V31 
WINDW3CW3) 
LINE3(W3(1 ) 


CAT.L 
CALL 
CALL 

CAT.L wiNuwj ; ,, 

CAT.L LTNF3(N3(1 ),W3(4) ,0 
CALL LINFR3(0.,-W3(4),u 
WC l=W3(l)+4.0 
WC2=W3f41-100. 

CALL SOFCTLCO) 


,0? 

O 


' 5 


* A - v . 


LftUU our CXI# l, v 3 . 

CALL CHARC'LTKF TO C»A*9 t r . 0) 
CALL LTNF3CWCi,wC2-50..-°- r 
CALL CHARC'WTNDW STZFC 
ACCEPT 31. A1 
FORMAT*! A1 ) 

IF(A1 ,NE.'Y')GOTO 32 f>. 

CALL LTNE3(WP1,W02-100 ? ;*5 *t 
CALL CHAR ( 'GT V* ME W2 0 •%„ ,U y 
CALL LINE3C*C1,WC2-150. ' 

ACCEPT *.(W2m,I = l,6) 

CALL CLRDFV (5) 

CALL VPQRT3CV1) 

CALL WINOW3CW13 
CALL MARGIN 
CALL WTNDWlCwn 
CALL VPOPT3CV23 
CAT.L RCIAOf 30. ,2) 

CALL ROTADf-30. ,1 ) 

CALL OBJECT 
CALL IDEM 
CALL DIMEN 
CALL VPOPT3CV31 
CALL WINDW7CW3) 

WF=*3(4)*.S 0 ) 

CALL LTNE3fW3(1 ),WF,0. ‘ 

CALL CHARf 'WANT PEREA T ' 1 
ACCE°T l ? , * S 
FORMAT C A1 5 
IFf ANS.N'F, 'Y'JCOTO 
CALL VPORTICVn 
CALL WINDW3 (WIT 
r* n t n t 

CALL CLRDEV ( 5 ) 


100 


100 



SUBROUTINE OBJECT 
COHMOfi/ACCEPT/CH,CC,X,' 
CALL LINE3C10. 0,900.0,1 
CALL. LINER1£X,0.0,0.0, 
CALL LINER 3(0. 0,0. 0,7,, 
CALL LINER3C-X.0. 0,0.0 
CALL LINER3(0.0,0.0,-Z 
CALL LTNER3(O.O,-¥»0.Q, 
CALL LINER3CX,0.0,0.0,' 
CALL LINER3C0. 0,0. 0,7,5 
CALL LINER3(0.0,Y, 0.0,5 
CALL LINER 3 { 0. 0,0.0, ”2 < 
CALL LINER3(0.0,-Y,0.0. 
FIN«. 025 

T=Y»FIN 

T GIVES FIN THICKNESS 
N*TFTX{Y/(2,*T)) 

CALL LINE3(iO.,900.,0. < 
DO 10 T*J,N 

CALL LINER3(0.0,-T,0.0, 
CALL LINER3CX, 0.0, 0.0,1 
CALL LTNFR3C0. 0,0. 0,7,1 
CALL LTNFR3C0.0,-T,0.0, 
CALL LINER 3 (0.0, 0.0, -2, 
CALL LTNER3(-X,Q.0,0.0 ( 
CONTINUE 
TT=902.-Y 

CALL LINF3(lO.0,TT,0.0i 
CALL LINERKX, 0.0, 0.0,1 
CALL LTNPR3C0. 0,0. 0,7,1 
RETURN 
end 


Y,Z 

0 0 . 0 , 

1) 

:‘5 

fj* 

1) 

1) 

*1) 

,1> 


SURR 
rov.M 
01 wE 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


OUTINE MARGIN 
nN/TRANS/Wl,W2,W3 
SION *'lC6),W2(6),tf3(6) 
LTNE3CW1 n).Hl£4),0..O) 
LINER 3 (Hl(2),0.«0.,t) 
LINER 3(0. ,-V?l(4) ,0. ,11 
LTNER3 ("*W1 (?),0.,0.,1) 
LTNER3(0.,W1£4),0.,1) 


SUBROUTINE TITLE 
CQWMON/TRANS/Wl ,*2,W3 
DIMENSION W1(6),W2C6) ,N3C6) 

DIMENSION Clf9),C2(7),C3m 
CALL SOFCTL(l) 

CALL CSIZEVC15) 

DATA Cl /'OPTIMIZATION METHODS F QR ENGINEERING DESIGN*. / 
DATA C2/'For . COMPACT HEAT - EXCHANGERS *.'/ 

DATA C3/'Staf f Advisor : Prof H.C. Agrawal*. */ 

CALL LTNE3CW1 (1) .Wi(4),0.,0) 

CALL LINFR3£W1C2),0.,0.,1) 

CALL LINER3(0.,W1 (4),0.,1) 

CALL LTNER3£-W1 (2) , 0. , 0. , 1 ) 

CALL LINER 3(0. ,W1(4},0.,1) 

MI=200. 

WN=100. 

WY-W 1(4 ) -200 . 

CALL LINE3CWW,WY,0, ,0) 

CALL CHARCC1 ,45) 

CALL LINE3CWW.WY-200..0. ,0) 

CALL CHAR £C2 , 31 ) a 

CALL LINE3(WW,WY-400.,0. ,0) 

CALL LINE3(WW,WY-400.,0. ,0) 

CALL CHAR(C3,34) 

RETURN 

END 



SUBROUTINE DTMEN 
CONMON/TRANS/W1-W2-W1 
COWMON/ACCFPT/CH,eC,X,Y,Z 
CQMMON/EXCCHR/FFF 
DIMENSION W1C6),W2C6),W3(6) 

CALL SOFCTL(t) 

CALL LINE3tW2Ct)+iO.,W?(4)-5O.,0., 
CALL CSIZEHHO) 

CALL CSIZEVCtS) 

CALL CSHEAf.2) 

CALL CHARC 'OPTIMIZED DIMENSIONS*.' 
CALL SOFCTL(O) 

Dl*tf2(2)*.5 

D2sW2(4)-200. 

CALL LINE3CD! ,D2,0.,0) 

CALL CHAR £ 'WIDTH :**.') 

CALL CHARF(X,7,2) 

CALL LTNE3fDt,D2-50.,0. i 0) 

CALL CHARC 'BREADTH :=*.*) 

CALL CHARFCY.7,2) 

CALL LINF3CD1 ,02-100. ,0, ,01 
CALL CHARC 'HEIGHT :=*.') 

CALL CHARFfZ.8.2) 

CALL LINFKD1 ,02-1*50. ,0.,0) 

CALL CHAR C 'AIRFLOW :=*.’) 

CALL CHARFfCH,^ l) 

CALL LTNF3fDl ,02-200. ,0., 01 
CALL CHAR ( *WAtER FLOW :=*.') 

CALL CHARFCCr.0,1 ) 

CALL LTNF3CD1 ,02-2*50. ,0. ,0) 

CALL CHAR C 'EFFECTIVENESS :=*.') 
CALL CAHRF C ERF ,6,4) 


CALL L 
CALL C 
CALL C 
CALL L 
CALL C 
CALL C 
CALL L 
CALL C 
CALL C 
CALL L 
CALL C 
CALL C 
RETURN 
ENn 


0 ) 

1 



RSSJIi'.r ,3F TRADE -OFF, STUDIES 


I VOLUME 
I 


METHOD 


IPENAL 



I PEN All 

EXPEN 

FEAS 


! IPENAL 

I * 

1 EXPEN 


I FEAS 


I IPENAL 


EXPEN 



| IPENAL 
EXPEN 


I FEAS 


IPENAL 


21.684 

21.548 

21.700 

15.173 

15.162 

15.400 

11.723 

11.716 

12.030 

9.592 


9.606 

9.770 


8.155 


8.126 


.550 


7.125 



I » 

XI 1 X2 


0.9431 


1.2870 | 0.8509 

1 ■ : S 1 ' • : ! 


1 1.2170 0.9540 

! 1.2106 0.8171 


I I 

| . 1,4648 j 0.7949 

I 1.1650 I 0.7943 

i f ■ - j i ; ! : * 


1-1.1341 I 0.6777 


I 1.1286 l 0.8039 

i i • i. ' • . ; 


I 1.0947 j 0.72,3 4, 



FEAS 

! IPENAL 

I " 

I EXPEN 

I ~ 

j FEAS 

I IPENAL 

I “ * 

I EXPEN 
i 


i 

I FEAS 


I 

! E 




0.7763 


1 I 

I 0.78421 

i t 


I 0.82511 

i i 


.3 


6 


! 0.85711 


l 0.85701 

i i 



I I 

I 0.88001 

i i 


I 0.88401 

i i 


! 0 .8920 ! 

i i 



I 1.1111 f 0.7389 

i • i • t 


1 0.89691 



7.103 

i . ' 

! 1.0946 

i _ _j 

1 0.7190 J 

| 0.9503 


























&PPEIOIX - B 


J PROGRAMS PACKAGE LISTINGS 






ft 

V 

ft 

V 


ft 

w 

ft 

V 


ft 

If 

C 

C 

ft 

L. 

c 

ft 

C 


w 

ft 

w 


32 


305 

300 

303 

304 

577 

572 

555 

571 

596 

591 

592 

593 


595 

590 

581 

582 

583 


570 


ms i55R5urpi”oiji5iI«ii’THi"iAiGi’fHrcr55NrliMFrdi"MiNlMOM 

MOUrni RANGE FDR MINIMUM LOCATED 
»3U£.r=2 CONSTRAINT IN EFFECT 

MQ (JIT's 3 M3 IMPROVEMENT POSSIBLE IN THE OBJECTIVE FUNCTION 
NREP=1 MO. OF ITERATIONS INADEQUATE TO ESI HATE THE RANGE 
NRE.P=<2 RANGE ESTIMATED 

MRBP=3 M3 IMPROVEMENT IN OBJECTIVE FUNCTION POSSIBLE 

rr,m = initial step size 

TMAX = M.AS STEP SIZE ALLLQiMED 

AL, AF = COORD. AND FUNCTION VALUES FOR QUAD. INTERPOLATION 
Z=I M IT IAL AMD FINAL VALUES ALONG A SEARCH DIRECTION 
KFEA,S=l UNCONSTRAINED OR EXTERIOR PENALTY FUNCTION PRJBLEM 
KF£A,S'=2 UNCONSTRAINED PROBLEM USING INTERIOR PENALTY FUNCTION 
XFEAS-3 CONSTRAINED PROBLEM USING FEASIBLE. DIRECTION METHOD 
METHOD CHOICE. OF MULTIVARIATE SEARCH TECHNIQUE 
METH30=1 UNIVARIATE SEARCH 
MET HOD =2 STEEPEST DESCENT 

MET HOD «3 CONJUGATE DIRECTIONS E POWELL'S METHOD ] 

METHOD *4 CONJUGATE GRADIENT C FLETCHER REEVES METHOD 3 

MET HOD® 5 VARIABLE METRIC l PATVDON FLETCHER POWElL METHOD ] 

I OUT — TYPE OF OUTPUT DESIRED 

NOTE. :» HIGHER ORDERS HAVE HIGHLY SUPPRESSED BUT 

IMPORTANT RESULTS 

I gum OUTPUT AT ALL STAGES 

TOUT =2 OUTPUT AFTER EACH ID MINIMIZATIONS 

l5ur=3 FINAL OUTPUT OF UNCONSTAINED MINIMIZATION 

I0UT®4 OUTPUT AT ALL STAGES CONSTRAINED MINIMIZATION 

r3Uf=5 FINAL OUTPUT OF &ACH CONSTRIANED MINIMIZATION 

I3UT=S ONLY FINAL OUTPUT OF CONSTRAINED MINIMIZATION 

COMMON/ M I N / A A L 
COMMON/CON VEC/G 

common/in:.ovc/go 

COMMOn/ALWAYS/TMAX, METHOD ,KFEAS f LIM.NC f NI. I OUT fSPERT 

COMnOM/pENFUN/R.UF 

COMMON/NSW/FF.ACCf IGRAD 

COM MO. n/ALARM/KALl JAW 

DIMENSION XU!*) f S(50),AF(3),AL{3),Z!2) 

DIMENSION YX!25),G(25) ,GO!25) 

DIMENSION DG!25,25) 

<NA = l 
AAL‘*TT 

imJUr.GT.2)G0T0 555 
WRIT S;!i,305)KAL, JAW 
WRITE* 1,3 99) 

WRIJ2C1,3)!X! J),J=l,Nfi) 

WRlfE.!l,353)ZU) 

WRlTSiUf 3D4HS(J) t J=lfNN) 

FORMAT! /3 5! '-*)/' ITER AT ION NO f 
FORMAT ('INITIAL VALUES OF X! J) ') 

FORMAT! 'INHAl VALUE OF THE OBJECTIVE FUNCTION ' /I PS 1 4 
FORMAT! 'INITIAL VALUES OF S! J) '/9(1PS1 4.6) ) 

IF! N CO 555, 555, 577 
WRITE! 1 ,5/2} 

FORMAT! f I NTIAL VALUES OF THE CONSTRAINTS') 

WRITE:! l,3)!G3!J},J = l,NC) 

H = ZU) 

IF ! N Ci) 5 7 9 , 5 7 9 , 5 7 1 

00 596 J=I,MC 
G! J ) =30! J) 

IF! XFSA 5-2)590, 591, 591 
CALL COnSTH! NPENAL j 
IF!.n?ENAL-1)590,591),592 
WRIT S'(1,593)KAL, JAW 

FORMAT! ** X! J) 1 INFEASIBLE IM',Ii,'tfl ITERATION OF ' 

1 ,12, * CYCLE ♦') 

WRITS! l,i)!&0(J),J=l , NC ) 

CALL SrARf!X,NN,NC,TT,DG,IOUT) 

CD MIO US. 

IF !XF£AS- 2)581, 581, 579 
IF! IOU T-l )582, 582,569 

FORMA^/'.'M PENAL' ,2X,'XNN' ,6X, 'STEP SIZE' ,8X, 'UNADGMETED' 


,12,' CYCLE NO ',12/35!'-')) 


>6} 


:rnN',i6x, ‘design vector') 


1 8 X , ' F U N 
'•n T O 5 59 
IF!lOUr-l)568,558,569 
jo rmi .<?£;? d 




kpe:=i 

ALU )= 3 . 

af a ) sz c i ) 

1 

KN = 1 
KX.= i 
KWH- 1 - 
K 21 =1 

fT=m " x,I4i>23,23?24 

?§iif t ,-f '*3 , 2 F ITERATIONS XOS THE MAX. SPECIFIED,- 
rDFN&r 1 ?* M3.0F JEXERATIONS INADEQUATE TO bOCATE THE MIN, 

£pf!u!*3* VIOLATED AFTER INITIAL DECEASE IN FUN^N. 

13? >> J ?iI^HLy I0LATfc:D RIGHT PR 3 « START.' 

SiU!Mi f 


JS(Jffl 51 4{5l4,SSl 


MC 


30 133 J , 

30 (JO *3.(13 

r* n C | A 

IteHs— , 

S 5 WlSi?l?' s, » 

m?A n 

30 l 3 1 J = 1 , m s< 
TX(J)=SCJ)+AAL*SCJ) 
:all> FiiN ? rc xxjF) 

CALL' M 331 ? Y(p) 

:ALL 33NSTH( NPENAL) 
|FCI3.0^-I J?56,535,557 
TFCMZjiSBl,bB4,5a5 
IFC£? iAWH3M87,584 
NRIX » 5 * 5 ) NPENAb.KNN, 

iinji ?f > r i ^ * .,?.» ? r r *> 2 « < 


£R| f fii A? * 5) NP£NAb,KNN,AAL,UF,F, UXl J) , J=l,NN) 
^J^^f Cl ^» 2 X),3ClPE14.6^X),5bciPEi4.6),/,&3X)) 

IE<^§AS-2)b4,6A,55 

rr*. 5 *rr ' ' 

AAL«TT . 

K 2 »*AZU 1 
30 TO IS 

SCHEME TO LOCATE THE CONSTRAINT. 

If c '[ :H055 r 66,S4 

30 o7 X=1 » Nr 
IF (3*C K) )6/, 58,63 

continue. ' 

3 = 3 C O 

Y| jj, 

30 S3 <=1,NC 
IFC SXiO )S3. 59 7 3 

\ J 

MC = * 

■’ON I I N J£ 

AAL = CSKMCH AftL - S (i (C )*HI)/(GOCMC)-G(MC3 3 

30 o J J.-l # M 

l*tl)=X .Ui+AAL*5C J) 

S.ALL- FJN ? |UX,F) 

CALu 403IFK F) 

wALb CONSTR^ENAb) 

If(L3Ur*l)^58558559 

H £‘ i s { iJUJii! i *’ 3 ' 75 


ivz c*r aT 7*1 f v,n 

lllil ikk ^ kklQ} 

3OT0 63 




104 


>3 104 
30CJJ*3.(j; 

7AU-F 
30 TO 122 
'KPSNAL = 2 

IFCK.FEAS-2) 82,82,83 

STS? SIZE CHANGE. CONSTRAINT 

OF 03U.EC,Tiye FUNCTION. 

rr=AAO-AAGO 

AAG=AAL'-TT 

F T = # 5 * T T 

AAL- AILM-TT 

KPE = 2 

KZl=KZitl 

GOTO 15 

I'F ( f *rt) SS, 55, 82 
4A03 s AA£i‘ 

XPESAL.*! 

1 ( 2 ) *F 

IFCZC2)-HJI4,14,15 

HsM2 : ) 

HI a A, A I i> 

IF.C 0)552, 552,533 

00 o5 j=i.nc 
30CJ0=3CJ) 

DETERMINATION AND UPDATING OF 
STEP SIZE SPECIFIED. 
fFCKZ -3)533, 530, 531 
AF(KZ)=ZC2) 

ALIKS) =AAL 
<Z=KS*i 
30 TO 532 
AFC 1) = AFC 2) 

AFC 2 ) = AFC 3) 

AFCi) = SU) 

AbC 1 ) = AuC 2) 

ALC2)=AuC3) 

A L I 3 ) = AAL 
CO NUN Jfi 

TF(<P£-I)49.49.48 
IFCKNn-33535, 537,537 
rr=AAL*CKNN+i) 

AAL=TT 
GOTO 50 
AAL=2'.*AAL< 

ST=AAL-AALO 

IF ( SF-IMAX) 51 ,51,52 

ST* MM 

AAL=AAL0+3T 

GOTO 51 

AAL = A’AL + Tr 

CONTINUE. 

K'NNs.KMNfl 
<K=KK>1 
<Z1=<Z1 5-1 
GOTO 15 

REVERSAL OF SEARCH DIRECTION. 

IF C KOI) IB, 19, 17 

IFCK.N-1)45,45,47 

GOT JC 23, 4 7, 2 0,47, 47 ), METHOD 

DO 19 J=1.NN 

51 J)=-SU) 

XN = K.N*1 
GOTO 15 

IF! K, £-3)533,533,534 

AFCKZ)=F 

ALIK.S)=AAL 

GOTO 535 

AFU) = AFI2) 

AF C 2 ) = AF C i ) 

AF C 3 ) = F 
ALC l ) = ADC 2) 

ALI2)=ALI 3) 

A L I 3 ) = A A L 
XL-ADiI 1 ) 

XH=AM3) 

MR5P=2 

HDUIF-1 

iF(rr-m)5oo,5oo,502 

rr=fTi 

GOTO 500 

IFCKNA-5)151,151,152 . 


VIOLATED AFTER INITIAL DECREASE 


AF AND AL. COMPARISON NITH MAX. 


DECREASE OF STEP SIZE. 


OA=*NA.+ l 
30T3 16 
M2Ui.r=2 
rr=x ri 

IFU:i) 153, 153,575 
>0 ?76 Jr!l , f}C 
3(J) =G3CJ) 

30 F 3 153 

MREP=4 

»OUI.Ptel 

if( rr^-rrn 50&,530, 503 

rr=rn 

30T3 500 

soui.r*.3 

IF ( MCI) 153,153,579 
DO 530 jsl ,«C 


3( J 3 s 3<3 ( J ) 
1(2) ‘ ' ‘ 


=2(1) 

NfRSP.sO ' 
fxZ ( 1 ) 

FORW(9(lPt'14.5)) 

ZUgtyM’* 013 561 

EXI1 no " ^ 

ifRfmi,85)MREP t NDUlr 

“ * “'If Sv * y'Gi.’."k » ^ 2 )( 9 j 


FORHATC 2X, ' MREP 
fRirs 
FORI AFC 
REF J.R*i 
END 


* \ £. n.0 intit ,«.<%» N0UIXV3X, X2,5X,I2) 
(l e 92)(G(J),J=i,ftC) 

" FRO'? RANGE AT EXI I ' /9 ( 1 PEI 4 . 6 ) ) 


ARE') 


44 

46 


78 

43 


79 

80 
81 


45 


7b 

82 

52 


83 

84 

85 


77 

8 

14 


9 

17 


23 
18 

24 


10 

11 


!H^25fiIi”i^§iIiIi^I^I^»n 5S=s?=aa=!!:aS!SS=s ' s,s ' =s=;sss=ss: = :s! 
:5Mpjiu , i5i"5F«lHisiuH‘'usISG'n54DRaTr" = fMTlDDnff5::5- J 

X -- 0SSI3.M VECTOR wuauraiu INTERPOLATION 


ORDEK OF THE DESIGN VECTOR 

&c.c accuracy desired a k 

S — . SEARCH DIRECTION VECTOR 

F-- FUNCTION VAGUE AI THE H2SIMUM 

GH 4AX SO. OF ITERATIONS WH 

AG COORDINATES OF THREE POLm 

At? FUNCTION VALUES AI THESE li&E POINTS 


■ I. ' L < 


J0'H;3S/PESFUS/R,UF 

§?8S/SKHK*‘ M '» eww ' , ' re “'M».'c,«i.»or.,?» I 


7‘Z 4T s**r 

cdmmom/mt n/aal 

“DMMOS/NEN/FF # ACC, I GRAD 
„OH*OM/AGARfl/RALi,$Afj 

21 SHIRS J^IhS(50),ALC3)aro) 


DIMENSION 22(25) 

IF ( LOUT-1 )44,44,45 




F0R4MC 1 RANGE RESJGTS FOP QUADRm" apRiyruaTny , 
*Rf|S;h,3)CAF(J), 0 = 1,3) AU “ *P R JXIMATI3,N > 

■J D f 1* f t 2 W U f 7 1 If « « -J \ 


sRirs:a,3)(AG(j),o = i;3) 
FDR 4 ATC 3(1 PE 14,5)) 


IF (3 Cl) 76,73,79 
ifft IT 2:(l t 4j) 

IF(Xf|as-2)80,8D,78 

4RIISC1 ,81 ) 

FORMAjfC/?*, ' p' , 4X, 'SIEP SIZE'-Iy 'liNfin'-MPH pm * t, 

1 I 2 X , DESIGN VECTOR') imU„.M£NTED ,7)C, 'FUNCTION 1 

continue 

<:i = i 
<C 2 = 1 
NUMs-1 
H=AF(2) 



JJ O J.- I f H * 

ZZ(J.) = X(J)+AAG*S(J) 

CAGGt FUMCTCZZ.F) 

CAGE- HOJIFYCF) 

IF(I0ur-l)76,76,77 
IF( VC!) 82,82,83 

NRIIECI, 52JNUV1, AAG,F, (ZZC K) , KsUtn 

£3RM.Ar(k,I2,2hP£f4.6 ) ,4CniPE{I:6),/,32X)) 

?FCLF£AS-2)84,84.82 

VRlTE'(l f 83)MU5!,AAG f UF,F,(2Z(I),r:i,» n n 

THE THREE POINTS. ' IIS ’ SMUFFLIW* OF COORDINATES OF 

IF ( A.3SI F-H ) , GE. ACC ) GOTO 7 
IF (NJH-GI 4)14, 15,15 
NJM=NUM+1 

IFC A*-(2j-F)9,55,10 
H = A F ( 2 ) 

IFC AAG-ALC 2) ) 1 7 , 7 , 1 8 

AF( 1 )=F 

&LC 1 ) = AA.G 

IF CLC 12-3) 2 3,25, 25 

LC1=1 

XC2=<C2+1 

30 TO 2 

AFC3)=F 

AG(3)=AAG 

IFC LC11 -3324,21,21 

XC2 = 1 

XC1=XC1+) 

30 TO 2 

H*S» {T 

IF(AAL-AG(2))11,7,12 

AF(3)=AF(2) 

AF(2)=F 
AG C 3 ) = AG ( 2 ) 

AG(2)=AAL 





ALC J.),rYCJ),(2Z(K),K=l,NM) 


computation of the mid-point too confirm 

SATISF rCATION F3 ACCURACY CRITERION 

00 / 2 K'Z 1 , M M 

ZZ(K)sX(K)tABSCALC3))*S(K) 

: a O'O' fonctczz.f) 

CALL MODIFY CFj 
YYC X) = F 

IFCf 3Ur-l)82,82,83 
IF(Om,?6,97 
NRlTSiC 1,52) A Li 
30 TO 83 
IF( (CFE.AS-2) 98 ,98,96 

SCLffAWCZf 

30 r 3 L p . 

AAL«ABS(AL«C3)) 

IFOft 2)-YY(l ) ) 76 , 76 , 77 

j =2 

SOTO 79 
n(3)=mn 
0 ?? 1 

DO 79 tC,= 1 ,Mn 
ZZ(K)=KIK)+AL( J)*SCK) 

SRU’Stl.n ulM.KAL 

FORMAT! '* ACCURACY NOI ACHIEVED TN ',13,' 

1 OF ',13,' tn l.D CYCLE *') 

00 7 3 J=l f MN 

SCCJ)=ZZU) 

F = YT C 3 ) 

RET JR H 
S'VO 


MINIMUM: AMD- 


ITERATIONS 


x’--"3i5i35"vicroR 

V N — ' ORDS.R OF THE DESIGN VECTOR 
MINIM -.-. CHOICE OF l.D SEARCH TECHNQUE 
MINI 4-1 GOLDEN SECTION 
nU*-2 QUADRATIC INERPOLATION 
rr INITIAL STEP SIZE 
ACC ACCURACY DESIRED 

F — OPTIH.UM VALUE OF OBJECTIVE FUNCTION 






CD <4 HON / I N CD VC /GO 
COMMOM/CONVEC/G 

COMMON/ALN AYS/TMAX, METHOD, KFEAS,LIM , NO, IQUT,NPERT 

COMMON /NEN/fF, ACC, IGRAD 

COM MON /A L ARM/ KAL, JAW 

DIM MS i:Off K C MN> , AL ( 3 ) , AF ( 3 3 , Z ( 2 ) 

DIMENSION SC50) ,G(25), 00(25) 

OK s i 

m*f* 

KALsO 

CALL< FUNCTCX., r ) 

5 ALL 1 MODIFY ( F) 

Z{| ) Sf 

ire*: i) 7 /, 77,7a 
DO ?? JsI.mC 

gou)=gu 5 

continue 

N=1 

kal=<al+i 

S(N)=1 . 

DO 11 J=1 ,Mm 
IFC J-SJ12, 11,12 

st y)?o. 

CO Vi I h J E. 


CALu RANGE! NN , N QUIT, NREP, TT,TT1 , X, S , AF , AL , XH, XL, Z) 
core (/I ,30,31) NREP 
~JVJ CiD. 21,21) NDUIT 
GOTO (9D0,D01)M1.NIM 
CALL. GO LO(X,S,NN, XL, XH,F) 

IJFJ 93J 

CALL' QUAD!X,S.NN,AL, AF,F) 

IFCOJ83, 80,81 
DO 3 2 J.«l , NC 
30 ( JO =30) 

Z!2)=F 

IF!«CTK-3)tf3, 84,84 
IF( A3SCZC 2 ) -Z( 1 ) ) - ACC)84 , 84 , 21 
MsV + l 
7 , C 1 ) =Z ( 2 ) 

IFCN-NN)22,22,25 
rFCM«-l)20. 20,32 
IF(NQUO-i)27,23,2D 
CONVERGENCE rule 
IF ( li ?E.R T-l )20,2 3,85 
CALL C0MVR3(X.S,Z.NN,NC,KTK,HK,F) 

GOTO C2/.36,1§7,20)KTK 
XALsCAL+i 
GO T 3 a 1 
PRINT OUTPUT 
F = Z ( 2 ) 

IF ( IOU T .31. 4 ) GOTO 7b 

F3rU'tC*pJWEL"S UNIVARIATE METHOD USED *) 
IFUOUr.Gr.SJffOTO 7b 
VRITSCI ,1 7 7 ) K A L . , 

FOR MAT! * N 0 OF 1.0 CYCLES ' ,13) 

IFU0UT.GT.3)GOr3 7b , 

NRI f E ( 1 , 2 1 1 ) !Xt J), J®1 ,NN) x 
FORMAT! 4 £ X ] := ' , RUPtl*.o) ) 

■ i n r 1 1 ■my? <* it % *f f \ ft. 


NRlIEtl.4) Z! 2 ) , ACC 


FORMAT! * [ 
F=Z ! 2 ) 

rr=f ri 
RE TURN 
END- 


,1PE14.6//'ACCURACY ACHIEVED ',F9.5) 


I55Z£iriIIIl£llEil5l5ii“iI^2===== ; 

x " - -”"3 5‘i I.5.V " vec tor ■' “ “ ' 

«« ORDER DF THE DESIGN VECTOR 
MININ' --CHOICE OF i.D SEARCH TECHNIQUE 
nn.^sl GOLDEN SECTION 
MINI.M*2 QUADRATIC INTERPOLATION 

rr initial step size 

ACC -- ACCURACY DESIRED urg 

FF percentage change in variable desired for FINITE 

DIFFERENCE; COMPUTATION OF THE GRADIENTS 
I GRAD* — CHOICE OF GRADIENT TECHNIQUE 
IGRAQcl USE PACKAGE ROUTINE 
IGRADC2 USE USER RDUTINE 

F — OPTIMUM VALUE OF OBJECTIVE FUNCTION , 

COM 4PM *)'£ 5$f EC/S'* '• * 

COMMON /IN CO VC /GO 

COMMOM/ALNAYS/TMAX, METHOD, KFEAS,LIM,MC,NI, IOUT, MPERT 
CDM 4PN/NEV/FF , ACC , I GRAD 
CDMMPN/ALARM/KAL, JA« 

DIMENSION X l 4 n ) , A L C 3 ) , A F C 3 ) , Z £ 2 ) 

DIMENSION S(53) ,3(25) ,G0(25) 

KrK'=i‘ 
m = rr 
kal=d 

CALL* FU.CTCX.F) 

CALL- MODIFYCF) 

ECU =F 

IFCMCO 7 7,77,78 
DO U J = 1UC 
GOCJ)=GU) 

CONI’ TM JE 
K'AL- K'A L + 1 

cor j (700.701 JIGRAD 
CALL GRADICX,S,SUM,NN,F) 

G o r 0 7 3 3 

CALu GRAQUCX,S,SUM) 

^ = 5^RTCSUMi) 

DO 3)3 J = l,JiN 
SC J) S5CJ.J/B 

CALL* RANGE! MN » N QUIT , NREP, IT, TT1 ,X , S, AF , AL , Xti,XL> Z) 

GOTO Cil, 12,13) NREP 
GOTO U2, 41. 41)MgUlT 
GOTO ODOfPDlJMINIM 
CALL* G3L0U,S,NS,XL,XH,F) 

GOTO UD3 

CALL* UUADCX,S,NN# AL,'AF , F) 

IFC 4:1)83,30,81 
DO 32 1 = 1, fiC 
GO ( J ) = G C J) 

U 2 )-F 

CONTINUE. 

I F C K DC.* 3)514,43.43 
IFCA3SCZC2)-ZC1))-ACC)43,43,42 
ZCU =Z(2) 

F=ZC2) 

GOTO ID 

CONVERGENCE RULE 

IFC .** ?£Rf*l ) 41 , 41 , 44 

CALL. COVVRGCX,S,£ t NN.NC,KTK,HK,F) 

G0r0ClD,43,14b,4UKTK 

KALsKAL+1 

GOT J IS 

PRINT’ OUTPUT 

r =Z ( 2 ) 

iFCiour.Gr.uGoro 7 & 

VRIT EC 1 ,2030 

FOR^TU STEEPEST DESCENT METHOD USED ') 

iFaour.Gr.5)Goro 73 
VRir SCI ,1 76 jKAL 

FORMAT ('NO OF 1.9 CYCLES U4) 

IF CiOUr.GT. 3)GOTQ 73 
«irai,3) C X CO ) * J S 1 » NN) 

NRlf2:(l,2J4) ZC 2) , ACC 


NRlf EH , 
FORMAT! i 
FOR MAT C ' 

rr=m 

R E T J R M 


•\ f i 


ACCURACY ACHIEVED %F9.5> 



F 2 = ( :IK *• F M 1$ ( M N ) "F M ) *2 
F3 = .5*FV|*CH.K-F)*2 
Zi=Fl*F2-F3 
TFUZm, 32,32 
50*=D. 

DO 3 5- 71=1, MM 
S'( JJSXCJJ -XO(J) 

SU^ s 'SU v f + S( J) *2 
3j*=saRr(suM) 

DO 35 ' 

5( J0=S(J,)/5UM 
IFOX«MM)40,4i,41 
MMX = 

DO 3 7 3*1. mi 

IF(X-MX)37,38,38 

DO 3? 1*1, HN 

500 ( 71,13 sSCOCJ+1,1) 

oom riwus. 

DO 42' J*4,«N 
soot W,X)=5(J) 

ISM 14 
fC «2 
' < 1=2 
fCAL«<’Ar.t 1 
'’OFO 25 

SUDFFLS r-iE DiREOriUM HATHIX 10 GET 
I = M M -1 

30 TO (33, 72,72), Kf 
IFC < r r- i -J) 14,12,12 
IF(Xri-^'U)&5,65,52 
IF(l-W:0 5l ,52,52 
4K = F 

zcu=f 

3orj r 33 , t 
DO 3 3 7 = 1 , MN 

50 D t *3 + 1 . j)=SOD(1, 1) 

DO 4 7 J = 1 , M M 
30 D(r, J)=S0DCI+1, J) 

DO 4 3 J = 1 , M M I 
DO 43 JX-1 , M M 
SOO tr + X, JX,i=SODCI + J+l , JK) 

<rr=<!|r+i 

33 ro 24 

OHEOiX; FDR C3MVERGESCE 

Ir (MP5RI-1)12,12, 73 

OALLi SON VRG ( X,S, * , NN,NC,KTK,HK,F) 

3Or0C?7,63,76,12)RrK 

DO 75 7=1 ,MS 

XOCJOsXl Jj 

KAL=XA5+1 

30 F 3 25 

DO 73 7=1, MM 

DO 73 1=1, MS 

SOD (I, 7) = 3, 

DO 7? 1=1, MM 
D3 7? 7=1 , MM 
IFtl-J) 79,80,79 
SODCt, J)=i. 
oomstiE 
<i=i 

< r r = i 
<:r< = l 
1 = 2 

< = 1 

DO 32 7=1, MM 
SODU, J) = S( J) 

30 r 3 24 

prist output 


F=Z ( 2) 

iFCiO 0 r. 3 r. 4 )Goro 173 

MRlf Etl ,58) 

FpR4 AI( l PjWELiL' 5 HEl'HQD t CDNDIR J 
tf ClDDT,Gr.5)GOrO 55 


SEAROH STARTED* 


USED* ') 



* * FC£rcHiR-iiiiis‘" 
F--"3i'ix5i'vicrOR""" s=sssses*sss»B»B*M*M«**«. 

™ •■:■ ORDER OF THE DESIGN VECTOR 
MINIM CHOICE OF |.D SEARCH TECHNI DOE 
MIMI.Msi GOLDEN SECTION ^ lWE 

J^?? A H^ INrE:Rt>OLA rxoN 
TT ** INITIAL STEP SIZE 
ACC -> ACCURACY- DESIRED 

SIr r -ri.E»?mE»IS e of l TS ,, iHllf ( if2»» P3R r ‘ virE 
iss&r^iPW^nffiiBi Rraii “ E 

I'GR AD«2 USE, USER ROUTINE 
f OPTIMUM VALUE OF OBJECTIVE FUNCTION 

****** * * * ' * - ’* '• - •« ■* • 

COM«S/I.NCOVC/G3 

com Mon/al^ats/tm ax .method, kfeas,l IM - Me# NI# I0UT, NPSRT 
COMMj3«/NSN/FF # ACC,IGRAD ' ' w ' " ,lJ “ i,,irwi 

C0MM3N/ALAR4/KAL, JA^ 


METHOD 


iftfi 


5),G0(25) 


DIMENSION XI MM) ,AF( 3 ),AL(i) 

DIMENSION AC 25 ),S 150 ) Sj(bO 

m*rr 

fCTK=1 
KA L s 0 

CALu F-MCiCX.F) 

CALu MODimF) 

IFCNC 077 , / 7 ,/g 

do n j=i f nc 

33C J)=3-( J) 

CONTI NO S. 

3K = F 

Z(U=F 

<al=<al+i 

GOT J I 7 DD, 7 DDIGRAD 
CALu GRAD! (X,S, SUM, NN,F) 

GOTO 703 

CALL GRA.DLJ ( X, S, SUM) 
isr-D 

3 = SJM 

SUM-SORT (SUM) 

DO 111 J=i,Ntf 
SOC JJ=S( J) 

5C J)=S( JO/SUM 

^OTO ( 51 s'S^sS } JIrep * ,TT,TTi # X,S,AF,AL,XH,XL', Z) 

mi fiftlfaiSiSitS 1 * 

CALL- 30 WHX,S,NN,XL,XH,F) 

DO 32 J= 1 ,MC 
GOCOJsGt J) 

CONTINUE 

UD-? 

IF I XTX-3 ) 5 3 , 92 r 02 
If CA3SU(2)-Z(I))-ACC)U,Il f 12 
rf(A3SCZ(2)-HK)-ACC)92,92 r 5i2 
MK=Z( 2 ) 

Z( 1 )=Z( 2 ) 

FsZl2) 

GOTO 513 
Z( 1 ) S Z( 2 ) 

F=Zt 2 ) 

IFUSr-,MM + l) 2 D, 20 , 5 l 3 

isr=rsm 

DETERMINATION OF CONJUGATE GRADIENT DIRECTION. 

GOTO ( 705 , 706 ) IGRAD 
CALL 3 RA 01 U,S, 5 UM,NN,F) 

GOTO 7 DA 

CALL GRAD(X,S,SJM) 

3 ETA-SUM /3 
3 —S JM : 

DO 45 J= 1 ,NN 
S(Jj=St J)t 8 ETA*S 0 ( J) 

S 3 (JJ=S£J) 

CO NT 1 21 US. 

SA = 0 .' 

DO 47 1 = 1 , NN 
SA=SA+S(I )»2 
DO 4 ? r=l # NN 
5 ( I ) = S( I ) /SORTC SA) 


30 TO 48 

::j;wsr3'£n:e criterion 

IF( M?SRl-l)5tl, 311,516 

CAUu CD. MVR3(X,S. 2, #H,|»C,KXK,HK,F) 

33 TO (5l3,lQ3,ll8,5inKTK 

■CAL*<AG>+1 

33 r 3 43 

PRIST DJTPUI 

F=2(2) 

IF(UUF,Gr t i)GOIO 76 

PDRWrFLSfCHER - POUELL METHOD (COSGRA) USED' } 

if( iDur.sf.sjGoro 73 

4RI ISI 1 , 1 76) KAL 

FORMAT! *W3 OF 1.0 CYCLES ',14) 

IFCiOUf ,GT« 3) GOTO 73 


<JRIf2;(l,205)ZC2),At_ 

FOR MAT C * t F 1 := ' , 1PE14.&/ 

rr»rn 

reform 

EM) 


ACCURACY ACHIEVED ',F3, 



A* 

c 

n 

w> 

ft 

U 

ft 

ft 

V 


103 


78 

79 
77 

124 


11 

12 

10 

700 

701 
703 


13 

65 

74 

16 

900 

901 
903 
81 
82 

80 

r* 

514 

515 
95 

121 

123 


90 

22 

705 


i £ OAVIDQffrFLETCH.ERrParfEUL' S_ 
x~—'~3isGLi~ vector'"”' 1 s=== 

MM - ORDER OF THE DESIGN VECTOR 

41 V 1.4; -- CHOICE OF 1,0 SEARCH TECHNIQUE 

4IMI4=1 GULDEN SECTION 

MIMX4s2 QUADRATIC INTERPOLATION 

rr — --.-initial step size 

ACC - — -ACCURACY DESIRED 

FF---;— -PERCENTAGE CHANGE IN VARIABLE' DESIRED FOR' FINITE 

difference; computation of the gradients 

I SR A.0— -CHOICE OF GRADIENT SUBROUTINE; 

IGRAOd USE PACKAGE SUBROUTINE 

IGRADcZ USE USER SUBROUTINE 

F— : OPTIMUM VALUE OF OBJECTIVE FUNCTION 

COMMO^foSvEa/Si'* * 

COM MON /IN CO VC /GO 

CDM MON /ALNAIS/T MAX, METHOD, KFEAS, LIM, NC ,MI , I OUT, NPSRT 

Sommom/min/aal 

C0MM3N/NE«/FF f ACC,IGRAD 
CDM43N/4LAR4/KAL, JA* 

COM40N/?EVFUN/R,F1 

DIMENSION <(NN) , AL(3), AF(33,Z(2) 

DIMENSION S(5j),H(25,25),FX(25),FT(25),G(25) 

SQ ( 25 } 


DIMENSION 

DIMENSION G3C25) ,G1 (25) 
START SEARCH IN STEEPDST 

r r 1 = rr 

kt< = i 

KAL=0 

CALL FJNCi'(XjrF) 

C ALu MODIFTCF) 
z ( 10 =f 

IFt-yCO 77,77,78 
DO n J = i , N c 
00 ( J ) = G1 ( J.) 

CONTINUE. 

RN = ? 

4 A L -■ I A L +■ 1 
DO 10 1=1, MN 
DO 10 I si , MN 
IF (I -0)1 2, 11, 12 

30^ I 3 ID* 

mm f 

GO TO (700,701) IGRAD 
CALL GiRADl(X,S,SUM,NN,F) 
GO TO 703 

CALL GRAOJCX, S, SUM) 
SUMsSQRT(SUM) 

DO 13 J=1 , MM 
FXC J)s-S( J.) 


DESCENT DIRECTION 


^(i,)=S(Jl 


)UI T , NREP , I'T, Til ,X,S,AF , AL , XH, XL, Z) 
t NREP 


. ;,tj =S( JO/SUM 
J ALL' RANGE! MM 

GO 13(514.16, . .. 

GO f 3(15,122, 12l)NQUIT 
GO T0(900,901)MIN1M 
CALu GOLO(X,S,NN,XL,XH,F) 

JJ 9 3 j. 

CALu JUAO(X,S,NM,AL,AF#F) 

IFCMCiJS®, 80,81 
DO 32 J=1 , Nw 
00(0) = G1(0) 

check* accuracy and restart criterion 

IFUrx-3)315.92,92 .. , 

IF(IOUr.LE.2)«RITE(l,95)(Z(I) 

-l O si ft IT f * if A - *“ “ “ ~ “ “* ‘ “ 


IF(iOUr.L£.2)«RITE(l,95)(Z(I),|slt2| ^ 

FORMAT! * VALUES OF ZC1 ) and Z(23 '/2C1PE1 4.6>) 
IF(A3S(Z(2)-Z(1))-ACC; 121, 121,22 

t tr m * ' nh nl 411 


IF( A.3S(Z(2) 

F=Z(2) 

HK = r 

Z(1)=Z(2) 
iFcnur.LE 
FORMAT!/ * 

GOTO 124 
s< — 2 i 2 } 

GOTO (705, 706) IGRAD 
CALL G'RAD1(X,S,SUM,NN,F) 

FOR^l 7 VARIABLE METRIC 


:)9i;92;i23 #22 


2)NRITfiCi,90) 

ITERATION IS RESTARTED 


FOR HEN SC JO »') 


11 ** 
11** 



06 

04 

4 


5frr*=?^' s ' slw 

oo ? § r.«i # »m 

33 3 3 1 = 1 , Mai 

S 5 IJ 4 .sSbLlM>HC I- J)*GCJ) 

F*CD = S5U* * 

ss 0 . 1 = 0 ,* 

corn's us. 
u=v 
03 5 2 ' J=U 

^2=Aa+soc J)»Gca) 

k 2 s kAi,y ^2 
00 61 I ® } # MM 
00 33 J = 1 . V M 

C i; r J.) +A2*S3( I3*SQC J) 

DO 5 ? X *1 , M S 

ft250+rrci)*3c j) 

00 >5' 1 = 1 , m N 
00 35 J = l, M s, 

P = t »- J >-FfU)*FXCJ)/A 2 

00 33 J = I * '•» ', 

FlCi) = -SC J.) 

SS'J'3 s O . 
oo 33 1 = 1, MN; 

03 *1 J = i , Hu 

S 5 aA = 3 SJ^rUI,J)*FifUJ 

S(I 3 =-S 5 JM 

socn = st u 

5 SUM =0 • 

**rj Ml HUS, 

Dfi'rs^i'i'An.DM 
OJ 32 1 = 1 * MM 

«■*.( j ) =? • i ( J) 

00 33 J = 1 , M 
5 S J^sSSUM + SC J )-2 

If(iiu 3 ) 67 ?l ^| 7 

S 5 iH = l. 

03 34 J’=l , M M 

IN 3Ur 880 " 

KAl,'=<Ab*l 
3 3 T 3 s >5 

CDS VSRjSJI CE CRITERION! 

IFCMPE.R 1-1 3122.1 2 2*516 

CALL* C 3 MVR 3 CX,S.Z,^M.WC,KfK,HK,F) 

PRIM I" OUTPUT 
p s t Z» C 2 J 

iK»Sf-fJ- 4,coro 75 

r 3R4AT( ‘OAVIDON - FLETCHER - POWELL 

rF(iour.Gr, 5 )Goro 76 

<j R I a ( 1 * 853 XAL 

n?«r(‘s 3 of i.o cycles M 3 ) 
IFCl 3 Ur.Gr ; 3 )GOTO 7 b 
/JRlfiU , 3 ) (XC J) ,J = 1 .NN) 

53 ?n?{ Uci> ,= •' 5(lPE,1 - ,,,) 

uitiy 1 F j 

rr=f ri 
M 3 = mOUU 
REIJ.RM 
SO 3 


OF VARIABLE METRIC DIRECTION. 


METHOD USED' 3 


= ' , 1 PE14 .6 ) 


this J 5 uii.irDirigicI 5 ii'fH£'’D£imiHis"uii 5 G"fii : "FlfiIir 
diffs.ks.mcs, method 

X -- QSSI. 3 H VECE.TOR 

MN — : ORDER of the design VECTOR 

F — INITIAL VALUE OF THE OBJECTIVE FUNCTION 

SUM *> SUM OF T«E SQUARES OF THE DERIVATIVES 

S SEARCH DIRECTION VECTOR 

£$«^Y-?S$M 5 *tt S !k*J!$SiHL lN VARIABLE DESIRED FDR THE 
CO MR J, I AT ION OF GRADIENTS 

■COMM; 3 M/AWWS/lMAX£MEfH 56 ^K^EAS^I*lM^c) 1 « JDUflH^SRr 

COMMiOM/NEN/FF.ACCjlGRAD 
DIMENSION X(NN),SC 5 Q),SS( 2 ) 

sscii=f nj 

SO.M’BJ , 

DO 10 Jfcl.NN 

IFUUD) 6 D, 61,69 
DaFf y ‘ 
pro hi 
D*FF*X(J) 

XC JO =Kt 03*0 
CALL- FUNCrCK.F) 

CALL MDDIFVCF) 

SS( 2 )-F 

X( J)sXC J .)-0 
sc imsscn-ssu) j/q 

SUM= 5 UMtSC JJ ^2 
IF (SJUil, 12,11 

SUM = I . 

CONTINUE. 

REFURv 
' J'*\i 3 


3.33 ROUTINE. GRADCX.S, SUM,NN,F,NC,DG) 

THIS ROUTINE DETERMINES THE DERIVATIVES USING THE FINITE 

DIFFERENCE. METHOD 

X — DESIGN VECETOR 

MM ->' ORDER OF THE DESIGN VECTOR 

F -- INITIAL VALUE OF THE OBJECTIVE FUNCTION 

SUM -- SUM OF THE SQUARES OF THE DERIVAriVES 

S' SEARCH DIRECTION VECTOR 

FF — ! PERCENTAGE. CHANGE IN VARIABLE DESIRED FOR THE 

COMPUTATION OF GRADIENTS 

DG — MATRIX OF CONSTRAINT GRADIENTS 

G, GO -- CONSTRAINT VECTORS 

NC NUMBER OF CONSTRAINTS 

. . . 

COMMON/CONVEC/G 

COM MON /I NCO VC /GO 

COMMON/NEN/FF,ACC t IGRAD 

DIMENSION X(MM),5C5Q),SS(2),DG(MC,MN) 

DIMENSION 30(25), 3(25) 

55(1 )=F 
5 J M = D . 

DO U 1 = 1 , MM 
IF(*( J>)53,61 ,63 
3 = FF 
GOTO »2 
3=FF *X( J) 

X(J)»X( J)*0 

CALu functcx, F) 

call MOOIFTCF) 

Si ( 2 ) =F 
X(J)=X(J)-Q 
S( J)=CSSC1)-SSC2))/Q 
■y 3 j i if ~ 1 . *i * 

OG(l,JO = (G.CXJ-G3(K)0/O 
SUM = SUM + SC J)**2 
IFCSJMm, 12,U 
5 J M = 1 . 

CON I IN US 

RETURN 

END 




|j|i5gr^ 


PJJJJJJJ 


:3*434/P54FU4/R,Uf 

:3^i,3,'i/::3'n/£C/,G 


£3M ID's / AO X A i’S/T MAX. METHOD r KFEAS , DIMM# Ml 

zd^osi/msa/ff.a:^, igraq 

DIME.'JSIJ'J 3(25} 


i fc *:«.sa.3)REruRN 
33rj(l,2# JXFEAS 


1 

(F(vljl3,13,14 

14 

33 5 J = 1 # 41 


T F ( 3.C J ) ) 5 » 5 , 6 

6 

FsFfR^GCjW 

5 

:j4r hue 

13 

4 £=.4 >41. 


I f C 4 £) 3 , 3 , 1 1 

11 

33 4 J=1,4E 

9 

F = F + R*3(NH + J) 4 * 


33 T 3 3 

2 

I F C 4 1 } 1 5 , 1 5 , 1 6 

16 

33 i 3- 1,41 


IF(3(J)}§#7 f 7 

8 

FsF-R/31 J) 

7 

:0*rEN35 

15 

4£=4 >41 


If('.£)3,3,12 

12 

33 1 3 J 5 1 f 4 E 

1 0 

F = c ’ + 3( 4jj>J)*2/ 

3 

R £ r 3 R .4 


£43 



17 

11 
1 J 

12 

14 

18 
10 

15 

16 


rc wpEmaEj" ”” = •==== s s;s:;i! 

J i ^ ^ k *■ 11 SwR S FOR V 1 OjLiA r 1 ON OF a N T ’’ 0 V S I S A I 'iTs 

NPiMAi,*! CONS fR 41 NTS SATISFIED -J'lilKAnri 

V?EN4.0 = 2^:,0NSrRAlNrs violated 

ACCURACY DESIRED 
F — J^aER OF CONSTRAINTS 

’“; l f V ^ A i ;|.-?r I ,5 A P!t? u ^ U ^E« Of VIOLATED 


riKr .. »f"rn3 » tu 7 r;% u cjuaiunfua vijyAi'sj 

NwJ'isi ve^ior aivTNG information on consirat nts violated 

x^^‘j^^jp^ 5S ^ A£Jfl55 ^ms;aM;Nc;Ni;iosi;5RiRT 

:O w 40N/CONVI3/N,,NCONSr 
wOM k *ON / NE N/FF# a5c # I SR AD 
OI^E.'iSlON G(25) , fl50NST(25l 

npeval=i 

IFt.<fSAS-l)16,l 5,17 
DO 13 J=1.NI 
IF(aCj)-,lS-loni ,12,12 
IFC A3SCGC J))-ACC)13,{ 3,10 
l iCD«srcM)=j ' ' 

SOTO 14 
NCD^Sl’I m) = J 
N P E * AL = 2 

IFC. < TEAS* 2)15, 15, 18 
N = M + l 

cdmhnue. 

30 To 1 3 
4 P £ V VO = 2 
Ri T j R * 


y.y if 



I n 

1 90 


93 

94 

91 

512 

100 

95 

83 


106 

107 

105 


99 

102 


101 

96 

97 


92 

96 


84 

85 

511 

80 


ES^ION/Edm VEC/G " 

:34 434/1.4:.3yC/G3 
Z'WM/WH/FF, a:z,igrao 

OIMS-MSHM 3(25) ,33(25) 

4X = Z (2) 

F-7A2) 

KU-U2) 

K'TK=2 

30 93 J=1,4M 
V( J)=X(J.) 

03 9# J=1 ,44 
U J)=1.01*X(J) 
calo- FUMcrtx.F) 

ZkUU :.34SfR( 4 IPE^KL) 

IFC «PS.f 40-1)80,80, 511 
Z2rfK-Z(2)-rfX 

IF(A35(Z2-K)-AC:)511 ,511,512 
If (4(2)-H!(.J99, 99,130 
33 95 3=1,44 
5(J)=AC J)-X(J) 

33 33 3 = 1 , 4 4 
X(3)=4(3) 

CALL F 3 4 C T ( X , F ) 

CALu OOIFHF) 

33 137 1=1,42 
30 ( J . ) — j C .J 3 
:39f 1.432. 

3(1 ) =3< 

33 tj m 

33 1)2 J = 1 ,44 

5C T)=XU)-A(J) 

4K=Z(2) 

Z( 1) =2(2) 

SJ.8-0. 

33 9 5- 3=1,44 
5d^=SimS( J)*2 
SU'l-SQRrCSUX) 

33 97 J=l,44 
5( 3)=S( JJ/SOM 
XTK = 3 
3313 80 

IF( A3S(Z(2)-HK)-A:C)511,511 ,98 
KTK-l 
F=Z(2) 

stilt i}” 

Z(( = rix: 

IF (M3!) 80,80,84 
33 65 3=1 , 4C - 
33 ( 3) = 3( J) 

33 r 3 80 

ktkm 

R'£ F J R 4 
Si43 





w 

f* 

w 

rt 

w 


302 

ft 

10 

15 

450 

451 
45 2 

453 

454 

455 


235 

76 

no 


35 

40 

13 


1 4 
116 

117 

144 


80 


ui|s‘rHi“iNT6Rl5R"’pilIXLr?’Fiiji5H55’rs5jiwI53i 

rj ? r si d Tie: constrained optimum 

X — . DS3IGH VECTOR 

N* -- ORDER OF THE DESIGN VECTOR 

VC — VJMjgR OF CONSTRAIN IS 

rr i vitial step size 

ACC — ACCURACY DESIRED 

*UU4 -- CHOICE OF 1 » D SEARCH TECHNIQUE 

MINIM=1 GOLDEN SECTION 
MlNiM=2 QUADRATIC INTERPOLATION 

METHOD — CHOICE OF MULTIVARIATE SEARCH TECHNIQUE 

M £ T H.3 D s 1 UNIVARIATE SEARCH 

METH0D = 2 STEEPEST DESCENT 

M£T,hDD = 3 CONJUGATE DIRECTIONS 

4£fHOD=4 CONJUGATE GRADIENT 

MET HOD = 5 VARIABLE METRIC METHOD 

<F£ 45=1 UNCONSTRAINED Ok EXTERIOR PENALTY FUNCTION PROBLEM 

<fe 4.5=2 unconstrained prolllm using interior pemauTt function 

<FEAS=3 CONSTRAINED PROBLEM USING FEASIBLE DIRECTION METHOD 
FF — PERCENTAGE CHANGE JN VARIABLE DESIRED FOR FINITE 
DIFFERENCE. COMPUTATION OF TtlF GRADIENTS 

R' -- penalty factor 

FI - -UN AUGMENTED OBJECTIVE FUNCTION 

NRE?=i TERMINATE. SEARCH DIRECTION DOESN'T IMPROVE. OBJECTIVE 

s-JNCTIOH OR VIOLATES A constraint 

FC -- FACTOR FOR INCREASING Ok decreasing "R* 




I m * m # • « m i 


CJM*lON/?E hfon/k,fi 
CD-* M ON/ ALN AY S/T MAX, METHOD, KFEAS # LIM,NC,NI, I OUT, NPEHT 
: J M 4 . 0 N / v s a / F F , A Z C , I G R A D 
CJ MmOn/ALARM/nAL, JAW 
CJMMO.N/COWEC/GI 
DI MENSI ON X ( 2 , 25 ) , CC ( 25) , DOC 25) 

DIMENSION V(2),RRC2),X(NN) 

DIME'S SI ON OGC25,25) ,3(25) ,Gl(25) 

J A N = l 

IF( i.oor.Gr,5)3oro 10 
VRl f S( 1 , 3D2)R» JAM 

FDRMAT(/4DC*- f )/.' R := *,tPE9.2, ' CYCLE NO ' . I 3/ 43 ( ' - ' ) ) 

PERFORM UNCONSTRAINED MINIMIZATIONS FOR TWO VALUES OF " R " 
NR = i 

GOTO (450,451,452,453,454)Mfc,THOD 
CALL UNIV(MI.SIM,X,rr,NN,F) 

GOTO 455 ' 

CALL' SrS£PCMINIM,X,TT,NN,n 
GOTO 455 

CALL' C,3NDIR(MINIM,X,TI»NN,F) 

GOTO 455 

CALL 1 CONGRA(MINIM,X,TT,NN,F) 

^3 fj 155 

CALL' DFPM(MXNIM,X f TT,NN,F) 

IFa.3UT.Lr.3.0R.I0Jr,GT.5)G0TC3 no 
«RirSi(l,235)(XCi),I = {,NN) 

<JRli’E:U,76)F 

FORM ATI * £ X ) ; , 9( IPE14.6) ) 

FORMAT! ' f F J := ',1PE16.8) 

CALL' FUNCTCX, F) 

CALS' MODIFY(F) 

s*i=F 

IFUOJT.Gr.SJGOrO 35 
VRl T EC 1 # 61 ) FI 
R R ( N R ) =u 
VI NKjsFi 
DO O J=1,NN 
VCNR,J)=Xi J) 

IF(‘i 4-2)13,14,14 
N R = N R f 1 

r=f:.*r 

JAft*= JA»tl 

TF(lOjr.GT.5)GOT3 15 
VRlYE(l,3U2)R,JAA 

** o r o 1 3 

BEGIN THE EXTRAPOLATION SCHEME 
IF( 4351 v ( 2) -VCl))-ACC)n6, 116,117 
F=VC2) 

GOTO 35 
MfAL : =l 

isSHSItSSliJJ 

0£M=SC-SD 

IF ( i.O 0 T. GT. 4 ) GOTO 70 
4 R I f S'C 1 , 4 6 2 ) 


4->? 
7 i ) 


lo 

8 2 
400 


401 

402 


403 

71 

2b 


17 


56 

53 

So 


8b 

85 

90 

93 

59 

91 

92 

87 

88 
89 

57 

405 


72 

456 

457 

458 
469 

460 

461 

111 


61 

73 

C 

37 

18 

19 


for m a t ( extrapolation results*) 

*A=( V(2)*$C-V(1)*SD)/DEn 
35=C/(1 )-V(2) 3/DEN 
3J i3 J = 1 

:"( J)=£n(2, J)*SC-N(1. J)4SD)/DES 
3J( J.) = UU, J)-*(2, J)5 /DeN 

tf’UDur.Sf. 4) GOTO 71 
NRIT 2(1 ,400) 

fjrmatc extrapolation scheme is r t i := kh + B 6 * R*.s v 
J predicted optimua is A A *) 

<VRI Ti(l .401 ) AA 
F 3 R M AT ( 1 A 4 ;= MPE15.8) 

<#8112(1,402) 

eormaxc ‘extrapolation .scheme is xaj=::(i)t3D(n*r,5V 
c^k'kSi.?" 1 "” aesla " ' ertot ls : - C15 ' 

rOR9Ar(9(lPElb,8)3 
IFCmL-1 325,25,26 

r=f:^r 

JA4=JAWtl 
A A 0 = A 4 

F M IS =4 A f B 8 * S Q RT ( R ) 

03 17 J=J,MM 

00 C J)=ZC( J)+D0( J)*S w Ri(K) 

CALL* FDmCI(DD,F) 

CALL MOOIFYCF) 

CALu C J4Sf«('JPE:iAu) 

I ft V?£¥4L-1)5o,S5,B5 
T F ( F 3- F 2)37, 57, 53 
0 J 33 3=1,11 
U 33=00 (J) 

ifu our.cr.4)Gor3 57 

*3RI i 2.1 1 , 8 a J R 

F JR 4.AT ( 1 EXTRAPOLATED VECTOR VALID tflTH NEW R 
3 0 T J 5 / 

CALu’ FJvCfCCC.F) 

CALu MODIFUF) 

CALL C O.lS TR( NPE N AL) 

IFC MPOAL-1 359,59,90 
IF(I3UT.Gr.4)GOrO 57 


, 1 P 2.9 . 2 3 


R I r 2C 1 f 9 3 > 
cor j 57 


FOR M AT ( 


REDICJEO VECTORS INFEASIBLE *'} 


IF(F3-P4}88,88,91 
03 92 J«1,NN 
XC J)=CC;U) 
liRlT SIC 1,87) 

F0R4A.K 'PREDICTED VECTOR CC VALID*) 

GOTO 5 7 
NR IT 2(1,89) 

F0R4ATC *EATRAPOLATED and PREDICTED OPTIMUM VECTORS INCREASE 
1 T HUE FUNCTION SO UNSUITABLE' 3 

iF(iour.Gr.5)Goro rz 

NRl i 2(1,405) 

F0RM4:r(*M£.*r starting POINT') 

4RI1S;U,433HX( J), J=1,NN) 

NRlTE(l,302)R»JAiV 

PERFORM UNCONSTRAINED MINIMIZATION FOR ME* VALJE. OF "R* 

GOTO (456, 457, 458, 459,460} METHOD 

CAL Li JUV(MIMIM,X, FT, IlN , F ) 

™ 3 T J i | 

CAL Li Sr£.£?(MlNIM,X,TT,Nh,F) 

GOT J 451 

CAL’y C0;iOIR(MINlM,X,TT,NN,F} 

GOTO 451 

CALL C0N3RAC M IN £M , X , IT , NN , F } 

7 3 F J 4 5 1 

call ofpmcumim,x. tt,nn,f) 

IF(10jr,LT.3.OR.IO<Jr,5T,5) GOTO 111 
<#Rir2Cl ,235)CX(I.),I = i,NN) 

VRlf2Cl,75)F 
CALL FUMCT(K,F) 

CALL 430IFXCF) 

iF(Uur.Gr.5)Goro 73 

/IRirSCl ,6UF1 . 

FORMA;T( f [ f ] maujHented := * , 1PE16.8) 

T*P C 4 4SC C 2 ) *F1 1 "13 11 ) 33 # M j 37 

-HSC<: ACCURACY CRlfi:Rl6N. IF MOT SATISFIED RE-COMPOTE CONSTANTS 

|F(A3S(Fl-FMliO-AC:)18,l8, 19 

NT AL=NTAL+1 

RR(U = RR(23 

RR( 2 ) =R 








i* 

r 


3 

450 

451 

452 

453 

454 

455 


7 

235 

UO 


Si 

1 

1 i 


L 

r* 

V* 

12 

5 


to 

4 

8 

2 


DESIRED FOR FINITE 


sub « 2^£^:L§*£iiil!!* N i«7x7tt7nn7^ fc~fJ 
^P-D^PyL^iFBP2il”F^’‘C0NSTRAiMiD”3FTl«!U' : i”ijiT^G*i : niF/ 

tAT^IOR PENALTY FtJMCIYDN TECHNIQUE 
X OiSIaN vector 

^ JRDS.R OF THE DESIGN VECTOR 
MINIM) -- CHOICE OF THE DESIGN VECTOR 
MIMI.Mst GJGDEN SECTION 
M I A i ,M = 2 QUADRATIC INTERPOLATION 
Tf INITIAL STEP SIZE 

NC. NUMBER OF CONSTRAINTS 

ACC ACCURACY DESIRED 
FF -- PERCENTAGE. CHANGE IN VARIABLE 
DIFFERENCE. COMPUTATION OF GRADIENTS 
IGRAD- -- CHOICE OF GRADIENT ROUTINE 
IGRA.Dcl (JSE PACKAGE ROUTINE 
IGRAQc2 USE USER ROUTINE 

METHOD — CHOICE OF MULTI VARIATE SEARCH TECHNIQUE 

METHODS! UNIVARIATE SEARCH 

METH0J=2 STEEPEST DESCENT 

M£ThQ0=3 CONJUGATE DIRECTIONS 

M£THD0=4 CJvJuCATE GRADIENT 

M£F H.3D-5 VARIABuE *£TkIC nEIHOD 

KFEASsi I.HCJND TRAINED OR EXTERIOR PENALTY FUNCTION PROBLEM 

kfsa5=2 unconstrained problem using interior penalty function 

KFEAS=3 CONSTRAINED PROBLEM USING FEASIBLE DIRECTION METHOD 

R -- PENALTY FACTOR 

FI -- UvAJGMSNTQ objective FUNCTION 

5-C — -‘'ACTOR FOR INCREASING OR DECREASING " R 

C ON 4 5 \ /P-C NFUN /R j FI * ” * * * V. V. v. , 

:0'*ON/At.«AYS/TNAX,METHOD,KFEAS,LIM r MC,NI ,I3UT, NPSRT 

:jMAO\/as*/FF,A:c , igrad 

commom/alarm/kal.ja* 

DIMENSION X(NN),CC(2) 

J=l ; J ANSI 

rFCiour. 3 r. 5 )Goro 3 
WRITE* l , 1D3R, JAW 

PERFORM UNCONSTRAINED MINIMIZATIONS FOR TWO VALUES OF R 
GOTO C 453,451, 452, 453, 454), METHOD 
CALL. UNIV(MIrNlH,X,rr,NN,F) 

GOTO 455 

CALL. SrE£P(MINIM,X,XT,NN,F) 

GOTO 455 

CALL C 3ND IRC MINIM, X, IT, NN,F) 

33 f '3 455 

CALL. C0N3RA(MINIM,X,TT,NN,F) 

'’□TO 455 

CALL OFPMC MINIM, X,TT.NN,F) 


JWUU Ur |r “ l H i ■'* A ~ , A , * * .I 

IFC l.OUT.Gf. 5 )GOiD UO 
»Rlf 2 (l, 7 )CXCI),I=l,NN 3 
WRlTEiCl J 235 )F 


, 9 C 1 PE14 . 6) ) 


(iurbii>< 

, 1PE14.6 J 


1 , 1P£1 6.8) 


FJRM4TC * C X 3 
FOR M.AT ( * I F 3 *s. 

CALL. FUNCTCXfF) 

CALL MOOIFYCF) 

:ccj.)=fi 

i?ciour.3r.5)Goro l 

WRlfECl.Si >CC(J) 

FORMAT! C F 3 uiatmtanted 
IFU-2)U ,12,12 
R = FC.*R 
JAW= JAWtl 
J = 2 

IFUOUr.GI.5)GOTO 3 
*M[fStl,lD)R,JAN 

^ J £ J j 

'H5C< ACCURACY CRITERION. IF NOT SATISFIED PERFORM UNCONSTRAINED 
MINIMIZATION FOR NEW VALUE OF R 
I F C A 33 C CC (1 ) -CC C 2 ) ) - ACC ) 4 f 4 , 5 
— C L )=’"'(2) 

R=FCTR; JAWSJAW+1 
I F ( 1 0 U T . G £ . 5 ) GO T 0 
WRIT S:C1,1D)R, JAW 
forma:t(40C'- s )/r 
GOTO 3 
F=CC.C2) 

iF(iojr.Gr.5)Goro 2 

WRIT5(1,7)(X( J), J=1,NN) 

#S5iSi‘ , [ , ? : } , Ss. J * lfS6..„;i0 gr 

1 /'ACCURACY ACHIEVED * ,1?F8. 

return 

S'ND 


3 


, 1PE9 . 2 , ' CYCLE: no ' , 12 ,/45 C ' *f )> 


constrained 

. 2 ) 


CrCLES ',13/ 





w 

C 


70 

701 

703 

111 

24 

U 

900 

901 
903 
40 

10 

22 


12 

13 

BO 


81 

83 

82 


79 

33 

31 

35 

705 


1^25^ ASC^l^XM £ x7fi7NN7w7Nc755.’J =- ==~=*^ 

f2i s 4rK i >iF 5 Pii|jijPofifii||B” 5 M”pn;5r35nr 

X -- OSS 1.3# VECTOR h DIhEa - T1DnS 
C-' 3RDER 3f THE DESIGN VECTOR 
41 '1 1 4 -- CHOICE OF l.D SEARCH T ?"4« r Tinr 
Ml MI H = 1 SODDEN SECTION Te.vdWI3UE 

0 ’5R?^l I ^ tNr ^ PDL ArioN 

feYH*u|fM E s i!a8{?2? Mtf ”“ e 

I li-f, JM.5W 

ACC -- ACCURACY DESIRED 

pp?E ? P|p5SiSrip?g%PGSiig5IS!fl 3ESIRSD FJ * ««« 

* R w-T OR 0? RUSH — Qi F F AC 1 0 RS 
MC *• MUM3ER OF CONSTRAINTS 

^OMoEr OF CONSTRAINTS VIOLATED 
NuONST -- VECTOR alVING IaFJ ONI CONSTRAINTS VIOLATED 
DG MATRIX Of CONSTRAIN T GRADIENTS ^ *“ J 

S', 30 -- CONSTRAINT VECTORS 




Awwf IGBUD 

iiil!il»EP i “« a '"“ = -" 


AF,AL,XH,XL,Z) 


CO 4 40 X / c J '4 v 1 3 / !i , N c 3 M S T 
CO4M0N/COWEC/G 
CD 4 4,3-4 / IN CO V C / GO 
CO^Ov/NEx/FF, ACC, IGRAD 
DIMENSION X ( M N ) “ e 
DIMENSION S ( 5 0 ) 

31 4 £, MS I ON A ( 2 5 
<0 = 1 

CALL f JvCT£X,F) 

7AU =f 

DO /O J=1 , NC 
C0( J.)=3( J) 
r r i = r t 

CO TO (703,791)1 SR AD 
-ALL- ORAD(X,S,SUM,MN,F,NC,DG) 

-jtO r J 7 3 3 

CALL* 3RADLKX, S, SUM) 

5U4*S0RT(SUH) 

32 111 0=1, <M 
Si J) =5C 0) /SUM 

CALL' RA.N3i.(NM.NDUir,NREP,TT,Tri,X,S, 

3-3 TO (10,11,12) 9ftEP f ' ' ' ' 

30 TO (900, 9D1)M£«IM 
CALL- 30LQ(X,S,MM,XL,XH,F) 

30 TO 9D3 

55^D^J?^5 s ^ j ^ AL ' AF#F) 

3-3(J.)=3(J) 

Z(2)=F 

ll\\^Yk\ U " Zil)}mKZZn22 fyi2,22 
‘*(17 “*.121 
F=Z(2) 

30 T 3 23 
F = Zt 2) 

30 TO (122.1 J. 122) JQUIT 
IF (<1.-2)83,81, 81 
F1=Z(2) 

<J = 2 
30 TO 82 
F 2 = ZU) 

IF(A3S(F2-F1)-ACC)122,122,83 

r l = r 2 
IX=M-l 

ill<; R '{3i4R I pR§GR»«5iNE OB ' >UIE “ SABW «»“•“ »*«*»* 

Ml = n +• M M 
M2=mM+1+M1 
M3 = n M 1 1 
DO 79 1=1, Ml 
DO 79 J=1,N2 
A ( I , J ) = D . 

DO 31 I = 1 , MM 
D3 31 1=1, MM 
IF(1- 0)31,33,31 
A(IA + UI, J,) = l . 

COM TIM US 
DO 35 1=1, MM 
3(IXH+J.)=2, 

30 TO ( 735,7 96), IGRAD 

CALL* 3RAD(X,S,SJM,MN,F,NC,DG) 

30 f 3 7 3 4 


> t « • # « * 



ZAlU' SRfcOJtX, S, SUM) 
33 3) J=1 , i X 
<=i:3Msrc j) 

3J 51 1=1 , 4M 
|CJ^n = 33(K # I) 

33 52 1=1 , MM 

:3Ml£MU£ 

33 55 J=1 ,IX 
K = M:,3sjSrCJ.) 

4( J, V3)srt(K} 

:3m,m. 

IX=IX+i 

A(rx,,U)=i. 

33 ?4 1*1, fil 


57 

IFU-J) 54,57,54 

\ci,mj)=t. 

54 

:3*m,vus 

33 71 J=l,42 

74 

:u)=o. 

«im = 3 

33 33 ) = 1,\>M 
4(U,J)=-$CJ) 

30 

5JM=§JM-S( J) 
3(lU = 5tH 

2 &Uu Li I 4 E A R ( A , 8 
73 7 3 J = 1 , M S> 

7 e 

5(J)=S(J)-1. 

SJ*=3. 

33 3? J*-l ,M<4 

59 

5u* = 5iM + SU)*2 
SJMsS^RTCSUM) 

I F (5 J '4)51,62, 61 

62 

5 JM = 1 . 

61 

33 5 3 J.= 1,MN 

60 

5CJ)=S(J)/SUM 
33T3 24 

122 

rfRirs(i,5)(xcj) 

F3RMATC 1 C X ) ; 
MRlISU .4)2(2) 
FORMjfciTl * [ F 3 S 

5 

4 

73 

REFORM 

2ND 


» 1PE14.6) 


V 

ft 


V 


101 

95 

60 

f» 

w 

102 

10 

97 

60 

94 

C 

87 

7 J 

89 

90 
72 


77 

91 
76 

92 
9i 


79 

80 


78 


85 

8 2 

83 

81 

74 

86 

56 

11 

12 

24 

44 


: -- »s5iui»fif§3||fiSffcgJ s M * M3 « SI3E DF Mi epK»tia»s 

P NIHBER 3F INEQUALITIES 

?.«« sf B |SE»B r «ftfss!SS a, “ R “ L * iUCK 

:0'<li54/AL«A^S/T4AX^ME!’HL)filKFEASlLI‘ir‘5 * ! NJ *f jfir *'i S? 5 1 ’ * * * * ’* * 

MRIi'2:(l ,95) 1,1 1,102 

Siiilfi ,s&'J SF*gjf £i GRA ' ,MI " G RES ‘ )MS ' 5 
[ 3 siA£S4tftii!Es 5 = i ! : r4r' I4 '‘“ os - ° F '“ lI * B1 ' SS 

FDR4. INITIAL BASIS 

00 n j=i , vi 

MBASISt J)=M;M-J 
:»( J)*C(V3+J) 

IFU3Ur-l )97, 97,8/ 
tfRII'Sd .63) 

FJROfC / 1 BASIS' , 5 X , ' CB ' , ttX , ' 6 ' , 60X , ' V ' ) 

3 w • J *® fc H # ^ i 

-§5iivj^ IJ:g8ASiStJ, ' C8(J),tt(J),UCJ,K),Kn, ’* 2) 

i T [, ^L a ,S ARE POSITIVE. IF N3T RESHUFFLE. BASIS. 

J J ( * J •* 1 f 1 1 

1 f ( ( J .) j 7 3 , 7 1 , 7 4 

IF ( 4.5'SC 5( J3 ) -i .£-10)89,89,90 
3 ( J J - 3 • 
sorj n 
90 n 1=1, M2 

au,<j=-au,k) 

3 C J J = -3 ( J J 

4 = 1 

03 /5 <=1,M2 

IF C 4C J, (0191.91, 77 

IF (At l, X) -1 . E-l 5 ) 91 , 91 ,92 

4 = 4 + 1 
CDNTOUE. 

IF ( 4-iM2)93,93,47 
3*A.U1,1) 

hc-i 

03 / 8 1=1, M2 
IFU.C J, (0)78,78,79 
I W-A.U., 0)80,78,78 

( Jlf K ) 

mc=i 

cSniNDE 


MBA5ISC J)=Hw 
:b(3)=:(h:. 


’*# »**' ’ '» W ■'£, W « TH » ^ 

03 35 1=1, M2 

StfjS?»fc K,/D 

33 61 <=1 ,M1 
IFC l-JO 82,81,82 
3=4(1, «C) 

33 83 1=1, M2 

A£<, £3=AC1,I)-D*AC J,I) 

B(0=B(0-D*B(J) 

COMI’IIMUE 

COMIOUE. 

33 35 J = 1 ,11 
IF ( 3 ( J ) ) B 7 , 86 , 8 5 

zsmaus 

S3 4=0 . 


3355 J=1 , Ml 

5 EL E.-ir COUdM T3 BE INCLUDED IN THE BASIS 
534=SJHt3( J) *CB( J) 

F = SJ.l 

33 12 <=1 , M2 
SL)H=0. 

33 11 J = 1 , U 
SU4=SUHtAC J,K)*:BU) 
zciU=su4-:;co 
:3 mii:mje. 

K.«=l 
M = 1 

33 43 1=1 , M2 
IFC4(iO-l. £-10) 44,44, 41 

« = M + 1 


4U 

41 

4/ 


30 

14 

li 


9b 

70 

71 
lb 


4b 

27 

28 
46 


31 

17 

1 6 

5 

9a 


61 

64 

63 


c 

5 


23 

18 

21 

20 

19 

55 

99 

75 

100 

f* 

43 

22 

50 

51 

103 

1 


13 


23 M £ L'M iJ £ 

IF(.*02)42,42,4 3 
y~i to 

33 13 J = 1 » M 2 

IF(iUJ)Ua3, 33 

5=?,ui <J,,u ’ 3,1 

vps J 

coin* us 

S£ti£.^l* R3/4 p BE REMOVED FROM THE BASIS. 

33 13 3=1 , Ml 
IFUCJ, 40)71 f 71 .95 

1 Ft AC J,4O-l.E-10) 71 ,71,70 
SCO =££3.) /.A CJ, MO 

a J X J 15 

cosxmo s. 

4 = 1 

33 d 7 3=1 , VI 
If CSC 33)45,28,23 

M'sMtl' 

continue. 

Iff 4-413 45,46,4 ? 

3=£ £3) 

4R=). 

33 15 1=1,41 
IFC5.( 3)315,31 ,31 
IF ( J-S.U) 316, 16,17 
3 = 2 1 J) 

MR=). 

2 ’3 4 I' 1.4 3 E. 

IFC<4-1 )5,S,6 
If U10r-1 3 93,98,5 
4 HI t' E {. 1 ,63) 

33 j4 3=1,41 

** '3 Nt f I s 4 U £ 

Jfftimi,6i)F, IZ(K).K=1,N2) 

FORM, ATC2X,'F = ',Flt).4,/2X, 'Z(K) = ' , 3 C 1 2F10 ,4 , / , 7 X) ) 

FORMATION OF THE. NErf SIMPLEX I ABLE B)f RSPLACInS R J4f IN f HE BASTS 
BY 1 HE. COULU4N. REPEAT THE CYCLE OF OPERATIONS. 

CBCMRlsCC 4C) 

Md45CSCMR)=4C 

3=Z(MO 

30 H K=1 , M2 

ZCK)=ZCX)-A(MR,X)*0/A(MR,MC) 

3= A t 4ft, HO 
30 19 3=1 , M2 
A(4R,))=A£4R, J)/0 
3(4R)=3(MR)/3> 

30 19 3=1, Ml 

IF().-MR)2l,19,2i 

3= AC j;, 40 

DO 23 <=1,M2 

AC J , O = A C ),, K) -D* AC MR, K) 

3(J)=3C3)-D*B(MR) 

CO Mf IN US 
SU 4 = 3 . 

33 55 1=1,11 

5J4=SU4t3C))*CBC J) 

?sS J.4 

IF CI3ur-l)99, 99,100 
4RITST1 ,53) 

33 / 5 3=1 , M 1 

\RliSCl,613'jBASIS(3),CB(J),B(J),CAC3fK),K = l,M2) 

CONTINUE 

MftlTSU ,6 3)F,CZCK),K=1,N2) 

<4 = 2 
con 24 
PRINT OUTPUT 
SJ«=3. 

DO 22 3=1, Ml 
SU4=SU4+BC))*CBC3) 
f = S J 4^ 

30 5 3 < = 1,M2 

ZCX)=3. 

30 51 3=1, Ml 
< = .M3A’SISC J.) 

Z C O =3 C J ) 

-nw 1 1 'H UE 

I F ( i OU T -.4)103,1 33,52 


RRIT5:C1,1)CZC J), J = 1.M3) 

FOR 4 ATC * C X J := ,9 C 1PE14.6) ) 



dRUIXl , i)F 

3 FJMVJri* l F 3 S=. MPfcl4,6) 

so rj vi 

4/ aRI JL S:(t ,4) 

4 FJFmi't' .43 S3L JTI 3 im POSSIBLE') 

52 oowdiue. 

REJ j.RM 
E4D 



nc 7 tt 7 dg 7 Iout) ’ • = - === = ===•=•= = = 

i'-ro 2 £ 0 S 3 r,?Et- r Si eEH G S vECtas 

I :> Slf fcfcKs s S88«« R0UIINES 

ff--- |,^f^ s !? E g^f£ge 

ixrrs*?S?l.^ 3 *p 5 jX?J 8 Sy"iHF , ' 5 MSfs 3 li II, * D rj * rm ™ 

DG HUTfUX 3F CONSTRAIN! GRADIENTS 


ooHov/coivtc/G* ' 

•'"'*» it 'll. 1 M ✓ if / £*• O Si **# if 




uomom/uodvc/go 

:o <5 fj.3 M /ms«/ff,a:g.igrad 

i &»**,««« 

UL ,J 1 s^ r J « s - ’sS^sfSS'MsKiMii'SSS!? 

fFCUjr.'ii.ifSuora 13 

tfUi SCl ,! 9) 

VRI f i c l , 1 3 ) 

** 1 x 20,1 ?mg(j), j=i,nc) 

*•5 = 1 

33 1% 3=1 , 4Z 
30( j)=G( J) 

a a l = r r 

33 3 K = 1 » iZ 
l F ( U * ) ) 1 ,2 , 2 
V = M + 1 
OO'iUtfJE 

TFC<-M:j7,7 f 10 

3=3(<) 

4 Z = < 

33 3 K=l, 50 

IF(GU)-5)J, 3,6 

3 = 3(0 

:5irr:sus. 

3 lfj R £/i!l, 7 lni 3 «AD RECTiON * , * 0RMftL t3 C3NS TRJUJir 43ST VI 
x^L 3RAD(X,S,S'JM,MN,F,l5C,DG) 

jJ I J / Ji 

CALL 3RA0U(X,S,SU.H) 

33 n j=i f f«M 

SC J3 =-£>3.C 
SUM =3. 

33 23 J=i,*N 

5U**5im5UJ*2 

su^ss^RfCsu*) 

IF£3JM)21,22,21 

53^=1, 

campus 

33 2 i J = 1 ,<5»5 
5(.TJ=S(J)/SUM 
iFti3ur.3£.2)Gora 15 
4Slf s:u ,35) 

FJRnATtOALUES JF S*J 

misu,ims(j), j=i,nn) 

33 12 J = 1 , *5 3 

fkU )=X( J)+AAL*S( J) 

CaLl, FJ'iCmx.F) 
iFU3ur.3s,2)5nro 93 
tRITEXl ,19) 

*Rll’E(l,17)(KX(J),J=l,NN) 

**£12(1,13) 


violated 


**£f£U,13) 

**U|(1,1/)(G(J),J = 1,NC) 

*A££2».(1,7 9) 

FDR^Ai’ ('VALUE OF 3(MC) and GO(MC)*) 
j 5 RllS(l f 17 )G(HC)#G 3 (MC) 

IFC 3>( *43 ) ~33( *10) ) 2 7 , 2 7 , 30 

30 il 1=1 f MM 

SfX(J.) = X£JJ+AAL01 E S(J) 

33 3 J= 1 , 5 V 

xcj)=mj) 

GOTO 23 
33 29 Jsl.SC 
33(J)=G( J) 

AAti3=AAL 


9 

2 7 
29 


14 

10 

94 

15 

19 

17 

18 

95 


dhfdx if all coNsrRAJNfa.fitfiL 'satisfied. if mot repeat me cstcle) 

IF(i.U:))32,UrH ■ 

DO Si 

KC J7=XXC J) 

GOTO 13 
AAl* = 2.*AAU' 
goto is 

IFCI 3 Ur. 35 .. 2 )GOrD 95 

mrs:ci t i5) 

FDR ^ATt 'FEASIBLE STARTING POINT') 
rfRtl SCI , 1 9 ) 

F 3 R 4 AIC ‘VALUES DF K ARE') 

dairs:a,i7)U(j) r j=i,NN) 

F 3 R* 5 Ar( 9 (iPE 14 . 5 )) 
tfftll SIC 1 / 1 8 ) 

F 3 JMATC VALUES 3 F 3 ') 

<Rif£;a,i 7 )(Gcj), j-i,nc) 

RETJ.RV 

END 



A 83^99 


fAS-' 






•V^OM 


D£-V 


